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A Mineral Survey for Piezo-Electric Materials 


By W. L. BOND 


those classes having centers of symmetry. 


propert ies. 


obtainable. 


the class in Groth’s Physikalische Kristallographie, as follows: (* 


ing classes of possible activity: 


*1 Assymmetric 


: ; Triclinic system 
2 Pinacoidal . 


Zeitschrift fur Kristalographie. 
? W. Voigt, Kristal physik. 
* Zeils f Physik 33, pg. 761 (1925). 
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ECAUSE of the increasing interest in piezoelectric materials in many 
branches of science an exhaustive study of the minerals was under- 
taken with the object of finding all the materials that could possibly be of 
use for piezo-electric elements. Much help was derived from existing data.! 
Considerations of symmetry show us that for a crystal to be piezo-elec- 
trically active it must belong to a crystal class that has no center of sym- 
metry (the Pentagonalicositetredral class of the cubic system, however, 
although it has no center of symmetry cannot be piezo active). 
makes twenty classes of possible piezo activity and twelve classes that could 


not possibly be active. About 90% of the crystals found in nature fall in 


Although the mineralogical data are incomplete in their assignment of 
minerals to definite classes in the seven systems, the existing data give a 
start in the choosing of minerals likely to have useful piezo-electric 


All available data were gone through to obtain the following list of min- 
erals classified by crystal structures. As many of the non-centric ones as 
were obtainable in the United States were tested by the method of Geibe 
and Scheibe* (resonance in a thermionic oscillator circuit). Whenever the 
authorities differed on the classification of a mineral it was so examined if 


In the mineral list, each mineral is numbered according to the number of 


'Dana—A System of Mineralogy, Ford—Dana’s Textbook of Mineralogy; Groth 
Chemische Kristalographie; Landélt Bérnstein—Tabellen; International Critical Tables; 





*3 Sphenoidal 
4 Domatic 
5 Prismatic 
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Monoclinic system 


*6 Bisphenoidal) 


* 


7 Pyramidal 


» Orthorhombic system 


8 Bipyramidal 


*9 Bisphenoidal 


*10 Pyramidal 


*11 Scalenohedral 
*12 Trapezohedral 
13 Bipyramidal 


*14 Ditetragonal Pyramidal 


Tetragonal system 


15 Ditetragonal Bipyramidal 


*16 Pyramidal 


17 Rhombohedral 

*18 Trapezohedral 

*19 Bipyramidal 

*20 Ditrigonal pyramidal 
21 Ditrigonal Scalenohedral | 


*22 Ditrigonal Bipyramidal 


*23 Pyramidal 


*24 Trapezohedral 

25 Bipyramidal 

*26 Dihexagonal Pyramidal 
27 Dihexagonal Bipyramidal 


» Rhombohedral system 


» Hexagonal system 


*28 Tetrahedral-Pentagonal-Dodecahedral | 
29 Pentagonal Icositetrahedral 
30 Dyakis-Dodecahedral 
*31 Hexakis-tetrahedral 
32 Hexakis Octahedral 


+ Cubic system 
| 


In addition to the above classification, the following list of minerals is 


annotated with the following symbols: 


nnd 


= active by test 
= inactive by test 
= unavailable or rare 

= mineral occurs only massive, amorphous or 


in other unsuitable form 


U = unstable 
C = electrically conducting 
? = class not absolutely certain 


S = crystal always very small 
{1 = mineral is always non-homogeneous 


CLASSIFIED LIST OF MINERALS 


Actinolite 
Adelite 
Aegirite 
Aenigmatite 
Aeschynite 
Alabandite 
Alamosite 
Albite 
Algondonite 
Allactite 
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Allanite 
Allemontite 
Allophane 
Almandite 
Altaite 
Aluminite 
Alunite 
Alunogen 
Amblygonite 
Amesite 
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Amosite 


Ampangabeite 


Amphibole 
Analcime 
Ancylite 
Andalusite 
Andesine 
Andorite 
Andradite 
Anemousite 
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\nglesite 
\nhydrite 
\nkerite 
\nnabergite 
\nnerodite 
\nomite 
\northite 
Anorthoclase 
Anthophyllite 
\ntigorite 
\ntlerite 
\patite 
\phrosiderite 
\phthitalite 
\pophyllite 
\ragonite 
\rdennite 
Ardunite 
\rfvedsonite 
\rgentite 
\rgentojarosite 
\rgyrodite 
\rrhenite 
Arseniosiderite 
\rsonolite 
Arsenophyrite 
Ascharite 
\strakanite 
Astrophyllite 
\tacamite 
Auerlite 
Augite 
Aurichalcite 
\utomolite 
\venturine 
Axinite 


Babingtonite 
Baddeleyite 
Baldaufite 
Barkevikite 
Barite 
Barytocalcite 
Bastnisite 
Baumhauerite 
Bauxite 
Beaverite 
Bechilite 
Beckelite 
Bementite 
Benitoite 
Beraunite 
Bertrandite 
Beryl 
Beryllonite 
Berzelianite 
Berzelite 
Betafite 
Bindheimite 
Binnite 
Biotite 
Bischofite 


MINERAL . 
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SURVEY FOR PIEZO-ELECTRIK 


Bismite 
Bismuthinite 
Bismutite 
Blédite 
Blomstrandine 
Boleite 
Boracite 
Borax 
Borickite 
Bornite 
Boulangerite 
Bournonite 
Braunite 
Breithauptite 
Britholite 
Brochantite 
Bromyrite 
Brookite 
Brucite 
Brushite 
Bunsenite 
Bytownite 


Cabrerite 
Cacoxnite 
Calamine 
Calaverite 
Calciothorite 
Calcite 
Caledonite 
Calomel 
Campylite 
Cancrinite 
Canfieldite 
Cannizzarite 
Carnallite 
Carnotite 
Carpholite 
Caryocerite 
Cassiterite 
Castorite 
Caswellite 
Catapleite 
Celestite 
Celsian 
Cenosite 
Cerargyrite 
Cerite 
Cerrusite 
Cervantite 
Chabazite 
Chalcanthite 
Chalcedony 
Chalcocite 
Chalcolamprite 
Chalcophyllite 
Chalcopyrite 
Chalcosiderite 
Chalcostibite 
Chamosite 
Chiastolite 
Childrenite 
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VATERIALS 


hillagite 
hloanthite 
hlorastrolite 
hlorite 
hloritoid 


‘hloropal 
hloraphoenicite 
hlorospinel 
hondrodite 
‘hromite 
‘hrysoberyll 
‘hrysolite 
‘innabar 
‘laudetite 
‘lausthalite 
‘leveite 
‘linoc hlor 
‘linoclasite 
linohedrite 
linohumite 
linozoisite 
obaltite 
ohenite 
‘olemanite 
ollinsite 
ollophanite 
oloradoite 
‘olumbite 
‘onnellite 
‘ookeite 
‘ordylite 
‘ornetite 
‘orundum 
‘orynite 
‘otunnite 
ovellite 
restmoreite 
ristobalite 
rocidolite 
rocoite 
ronstedtite 
rookesite 
ryolite 
‘ryolithionite 
uprite 
“upros¢ heelite 
‘yanite 
‘yrtolite 
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Dahllite 
Danburite 
Datolite 
Dawsonite 
Dechenite 
Delessite 
Dellafosite 
Delorenzite 
Delvauxite 
Demantoid 
Deschloizite 
Desmine 
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Deweylite 
Diamond 
Diaphorite 
Diaspore 
Diopside 
Dioptase 
Dixenite 
Dolemite 
Domeykite 
Douglasite 
Dufrenite 
Dufreneysite 
Dumortierite 
Dysanalyte 
Dyscrasite 


Edingtonite 
Eleonorite 
Ellsworthite 
Elpidite 
Embolite 
Emerald 
Emmonsite 
Emplectite 
Enargite 
Enstatite 
Eosphorite 
Epidesmine 
Epididymite 
Epidote 
Epistilbite 
Epistolite 
Epsomite 
Erikite 
Erythrite 


Euclase 
Euchroite 
Eucolite 
Eucairite 
Eudialyte 
Eudidylite 
Eulytite 
Euxenite 


Fairfieldite 
Fassaite 
Faujasite 
Fayalite 
Ferberite 
Fergusonite 
Ferrierite 
Florencite 
Fluocerite 
Fluorite 
Forsterite 
Forshagite 
Fouquerite 
Fowlerite 
Francolite 
Franklinite 
Freibergite 


Erythrosiderite 
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Freyalite 
Frieseite 
Fritzscheite 
Fuchsite 


Gadolinite 
Gageite 
Gahnite 
Galena 
Ganomalite 
Garnet 
Gastaldite 
Gay-Lussite 
Gedrite 
Gehlenite 
Germantite 
Gersdorffite 
Geyserite 
Gilsonite 
Gismondite 
Glaserite 
Glauberite 
Glaucodot 
Glauconite 
Glaucophane 
Gmelinite 
Goethite 
Goslarite 
Graphite 
Greenockite 
Griffithite 
Grossularite 
Guanajuatite 
Gummite 
Gymnite 
Gypsum 


Hackmanite 
Haidingerite 
Halite 
Halloysite 
Hambergite 
Hancockite 
Hanksite 
Hardystonite 
Harmotone 
Hatchettolite 
Hauerite 
Hausmannite 
Hauynite 
Hedenbergite 
Hedyphane 
Heintzite 
Hellandite 
Heloite 
Helvite 
Hematite 
Hercynite 
Herderite 


Herrengrundite 


Hessite 
Hetaerolite 
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3221 
*281 


*111 
*31I 
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Heulandite 
Hielmite 
Hieratite 
Hillebrandite 
Hiortdahlite 
Hisingerite 


Hodgkinsonite 


Hoeferite 
Hokutolite 
Holmquistite 
Hopeite 
Howlite 
Huebnerite 
Humite 
Hussakite 
Hutchinsonite 
Hyalophane 


Hydroboracite 
Hydromagnesite 


Hydrozincite 
Hypersthene 


Ilmenite 
Ilminerutile 
Ilsemannite 
Ilvaite 
Inesite 
Todembolite 
Iodobromite 
Todyrite 
Tolite 


Jadeite 
Jamesonite 
Jarosite 
Jeffersonite 
Jenkinsite 
Jezekite 
Johnstrupite 
Jordanite 
Joseite 


Kainite 
Kalinite 
Kaolinite 
Kasolite 
Kelihauite 
Kentrolite 
Kermesite 
Kieserite 
Klaprotholite 


Klebelsbergite 


Knopite 


Kobaltmanganerz 


Koenenite 
Koppite 
Kornerupine 
Krennerite 
Kroehnkite 
Kunzite 
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Labradorite 
Langbanite 
Langbeinite 
Langite 
Lanthanite 
Lapis-lazuli 
Laumontite 
Laurionite 
Laurite 
Lautarite 
Lavenite 
Lawsonite 
Lazulite 
Lazurite 
Leadhillite 
Lehnerite 
Lehrbachite 
Leonite 
Lepidolite 
Lepidomelane 
Leucite 
Leucophanite 
Leucophoenicite 
Libethenite 
Limonite 
Linarite 
Linnaeite 
Licroconite 
Liskeardite 
Lithiophilite 
Loeweite 
Loellingite 
Loparite 
Lorandite 
Loranskite 
Ludlamite 
Ludwigite 


Magnesite 
Magnetite 
Magnetoplumbite 
Malachite 
Malacon 
Mallardite 
Manganhedenbergite 
Manganite 
Manganophyllite 
Manganosite 
Manganotantalite 
Marcasite 
Margarite 
Margarosanite 
Margasite 
Marialite 
Marignacite 
Marmolite 
Marshite 

Martite 
Mascagnite 
Matlockite 
Maucherite 
Meionite 
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Melanite 
Melanocerite 
Melanophlogite 
Melanterite 
Melilite 
Meliphanite 
Mellite 
Mendozite 
Menilite 
Merwinite 
Mesolite 
Metacinnabarite 
Meta Torbernite 
Metavoltine 
Miargyrite 
Microcline 
Microlite 
Microperthite 
Microsommite 
Miersite 
Milarite 
Millerite 
Mimetene 
Mimetite 
Minium 
Mirabilite 
Mizzonite 
Molybdenite 
Molybdite 
Monazite 
Monticellite 
Montmorillonite 
Montroydit 
Morensonite 
Morganite 
Mosandrite 
Mossite 
Mottramite 
Muellerite 
Muscovite 
Muthmannite 


Nadorite 
Nagyagite 
Natrolite 
Natron 
Naumannite 
Nemalite 
Neotantalite 
Neotocite 
Nephelite 
Nephrite 
Neptunite 
Nesquehonite 
Niccolite 
Nickolsonite 
Nickelbluete 
Nickeleisen 
Niter 
Nocerite 
Northrupite 
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VATERIALS 


Noselite 
Nowmeite 


Ochrolite 
Octahedrite 
Okenite 
Oligoclase 
Olivenite 
Olivine 
Omphacite 
Onofrite 
Opal 
Orpiment 
Orthoclase 
Osmiridium 
Otavite 
Ottrelite 


Pachnolite 
Pandermite 
Paragonite 
Parahoepite 
Paralaurionite 
Paratakamite 
Paravavxite 
Pargasite 
Parisite 
Patronite 
Pearceite 
Pectolite 
Penninite 
Pentlandite 
Percylite 
Periclase 
Peristerite 
Perovskite 
Perthite 
Petalite 
Petzite 
Pharmacolite 


Pharmacosiderite 


Phenacite 
Phillipsite 
Phlogopite 
Phosgenite 
Phosphoferrite 
Phosphophyllite 
Phosphosiderite 
Phosphuranylite 
Pickeringite 
Picotite 
Picromerite 
Piedmontite 
Pinakiolite 
Pinguite 

Pinite 

Pinnoite 
Pirrsonite 
Pisolite 
Pitchblende 
Plagionite 
Plattnerite 
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Pleonast 
Plumbojarosite 
Polianite 
Pollucite 
Polybasite 
Polycrase 
Polydymite 
Polyhalite 
Polymignite 
Powellite 
Prehnite 
Priorite 
Prismatine 
Probertite 
Prochlorite 
Proustite 
Pseudobrookite 
Pseudomalichite 
Psilomelane 
Psittacinite 
Ptilolite 
Pucherite 
Pumpellyite 
Pyroargyrite 
Pyrite 
Pyroaurite 
Pyrochlore 
Pyrochroite 
Pyrolusite 
Pyromorphite 
Pyrope 
Phyrophanite 
Pyrophyllite 
Pyropissite 
Pyrosmallite 
Pyrostilpnite 
Pyroxene 
Pyroxmangite 
Pyrrhotite 


Quartz 
Quenselite 
Quercyite 
Quisqweite 


Ralstonite 
Rammelsbergite 
Raspite 
Realgar 
Rhabdophanite 
Rhodochrosite 
Rhodolite 
Rhodonite 
Rhomite 
Richterite 
Richardite 
Riebeckite 
Rinkite 
Rinneit 
Ripidolite 
Risorite 
Riversideite 
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Romeite 
Roscoelite 
Rosenbushite 
Rowlandite 
Ruby 
Rumpfite 
Rutherfordine 
Rutile 


Safflorite 
Sal-ammoniac 
Salite 
Samarskite 
Sanidine 
Sapphirine 
Sarcolite 
Sartorite 
Sassolite 
Scheelite 
Schefierite 
Schirmerite 
Schizolite 
Schorlomite 
Schreibersite 
Schrockingerite 
Schrotterite 
Schwartzengergite 
Schwetzite 
Scolecite 
Scorodite 
Semseyite 
Senarmontite 
Sepiolite 
Serpentine 
Serpierite 
Shortite 
Siderite 
Sillimanite 
Sipylite 
Skemmatite 
Skutterudite 
Smaltite 
Smithsonite 
Sodalite 
Sodaniter 
Spencerite 
Spessartite 
Sphalerite 
Spinel 
Spodumene 
Spurrite 
Staffelite 
Stannite 
Staurolite 
Steenstrupine 
Stephanite 
Sternbergite 
Stibiconite 
Stibiotantalite 
Stibnite 
Stichtite 
Stilbite 
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Stilpnosiderite 
Stolzite 
Strengite 
Stromeyerite 
Strontianite 
Struvite 
Sulfoborite 
Sulfur 
Sulvanite 
Sussexite 
Svanbergite 
Sychnodymite 
Sylvanite 
Sylvite 
Symplesite 
Syndalphite 
Synganite 


Tachyaphaltite 
Tachyhydrite 
Talc 
Tantalite 
Tapiolite 
Tarbuttite 
Tasmanite 
Teallite 
Tengerite 
Tennantite 
Tenorite 
Tephroite 
Tetradymite 
Tetrahedrite 
Thalenite 
Thaumasite 
Thenardite 
Thermonatrite 
Thomsenolite 
Thomsonite 
Thorianite 
Thorite 
Thortveitite 
Thuringite 
Tiemannite 
Tiger-eye 
Tilasite 
Titanite 
Titanmagneteisen 
Topaz 
Topazolite 
Torbernite 
Tourmaline 
Trechmannite 
Tremolite 
Tridymite 
Trimerite 
Triphylite 
Triplite 
Triploidite 
Tritomite 
Troegerite 
Troilite 
Trona 
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Troostite 17 Variscite 8?S Wollastonite 5 
I'scheffkinite M,H  Vauxite I Wulfenite *10 
I'schermigite 30 ~=Vermiculite I Wurtzite *20A 
lungstenite M __ Vesuvianite 15 
Tungstite 8 Villiaumite a . 
Turgite I Vivianite 5 re se J 
Turquois 2  Volborthite I Xanthos a ite = [ 
Tychite 32 Voltaite oor ong = 
Tyrolite M __sVonsenite anaes 15 
[ysonite 27 
Rae Wad MY ttrialite M 
Ulexite M : -e “erite 
see ae Wagnerite 5 ttrocerite M 
een ee Walpurgite 7 Laie eso 32 
Uraninite 32 Warwickite I seaecune : 
Uranocricite 8? Wavellite i I —— 
Uranophone M Wernerite 101 
Uranopilite M W hewellite »  Zeratite M 
Uranosphaerite M_ Whitneyite MI = Zaophyllite I 
Uranospinite 8? Wiikite I Zeunerite 15 
Uranothallite 8 Wilkeite I Zincite *1971 
Uranothorite M__ Willemite 17 Zinkenite 8 
Uranotile 2 Wilsonite I Zinwaldite 5 
Utahite ?S  Witherite 8 Zircon 15 
, “ Wittichenite 8  Zirkelite 32 
— °* — Woehlerite 5  Zoisite 8 
Valentinite 8  Wolfachite 8  Zorgite M 
Vanadinite 251 Wolframite 5 Zunyite *31AS 


Of the 830 minerals listed 70 belong to classes that allow piezo-activity 
but only 17 are found to be active by the Giebe and Scheibe test. (Our test 
of Iodyrite was negative but Greenwood and Tomboulian‘ found it to be 
active; on the other hand, we found Scolecite to be active while they report 
it inactive.) It may be that others of the remaining 56 classes have such 
small piezo-electric constants as to be undetectable. Others may be in- 
correctly classified as to symmetry. 

Of these active materials, quartz is the most important. Because of its 
excellent mechanical properties (stability, etc.) as well as for its relative 
cheapness it seems destined to remain one of the most important piezo ma- 
terials. 

Tourmaline is also important because of the high magnitude of its elastic 
moduli in certain directions; however, it cannot be obtained in large pieces 
of satisfactory homogeneity. 

Sphalerite is very difficult to handle because of its many cleavage planes, 
and appears to give little promise of becoming practically useful. Its activ- 
ity is quite marked. 

Homogeneous crystals of calamine appear to be very rare, so that work- 
able crystals large enough for ordinary piezo-electric application are un- 
obtainable. Most of the material occurs massive. 


*On Piezo Electricity--Greenwood and Tomboulian—Zeits. f. Krist. Jan. 1932. 
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Epsomite gives a marked response but the crystals are generally small and 


they do not weather well. There is some possibility, however, that they 
can be made artificially. 

Boracite gives a marked response, but boracite alters slowly. Its imper 
manence may bar it for some uses. 

Stibiotantalite occurs only in thin scales, and the necessary cuts must be 
made in the most wasteful way. Twinning is prevalent and the composi- 
tion varies widely. 

Scolecite occurs only as small crystals a few millimeters in diameter and a 
centimeter or so in length, uniformly twinned. 

lodyrite has been found to be active by other investigators. It is 
electrically conductive, very soft and not very common. 

Struvite is soft, unstable, and occurs only in small crystals. 

Zunyite occurs only in minute crystals. 

Langbeinite slowly changes its crystal structure. It may be made ar- 
tificially so may be of some use if it can be kept from alteration. 

Leucophanite and Meliphanite are related minerals. Neither seems to 
occur in good (i.e., homogeneous and untwinned) crystals of usable size. 

Wurtzite does not appear very active but good crystals were not 
obtainable. 

Tiemannite crystals were also unobtainable, but fragments of massive 
tiema nite responded. Crystals might respond more energetically if they 
were obtainable, but minerals that are too difficult to get would not be of 
practical use. 

Epistilbite occurs only in small specimens, uniformly twinned. 

The mineral clinohedrite is strongly active but crystals are very rare. 

Cronstedtite and Edingtonite are very weakly active. Crystals of these 
are very rare. 
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The Fundamental Equations of Electron Motion 
(Dynamics of High Speed Particles) 


By L. A. MacColl 


I. INTRODUCTION 


In work relating to the motion of electrons and other particles it is fairly 
common to assume that the particles obey the laws of Newtonian dynamics. 
That is, briefly, it is assumed that the rectangular coordinates (x, y, 2) of 
the particle under consideration satisfy the differential equations 


ma = X, my = Y, mi = Z, 


where m is the mass of the particle (assumed constant), Y, )’, and Z are the 
components of the applied force, and the dots indicate differentiation with 
respect to the time ¢. 

However, it is well recognized now that the above equations are not 
strictly correct, and that they merely represent an approximation which is 
adequate when the speed of the particle is sufficiently small compared with 
the speed of light. The system of dynamics based upon the correct equa- 
tions! (which will be exhibited presently) is commonly called relativistic 
dynamics, not because any knowledge of the theory of relativity is essential 
to its understanding and use’, but because it is in agreement with the theory 
of relativity (which Newtonian dynamics is not), because it was first de- 
veloped in connection with work on the theory of relativity, and because 
even yet virtually all of the expositions of the subject are to be found in 
books and papers dealing primarily with the theory of relativity. 

Just where the dividing line should be set between cases in which New- 
tonian dynamics is an adequate approximation and cases in which it is 
necessary to use relativistic dynamics is, of course, a rather vague question 
which cannot be answered simply and definitely. We may note, however, 

'It is not the purpose of this article to discuss questions of fundamental physics, or 
the physical validity of any particular equations. For purposes of discussion, we assume 
outright that relativistic dynamics is at least more nearly correct than is Newtonian 
dynamics. 

2 The theory of relativity can be described briefly as a theory of the relations between 
the descriptions of phenomena in terms of different systems of reference. We shall not 
be concerned with this theory, because we shall be employing the same reference system 
throughout most of our discussion. In the final section of the paper we shall consider 
purely geometrical transformations of the coordinate system. These transformations, 
however, involve nothing that is really characteristic of the theory of relativity in the 
usual sense. 
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that according to relativistic dynamics the mass of a five thousand volt 


electron is about one per cent greater than the mass of an electron at rest 


From this we can infer that, while Newtonian dynamics may be adequ 
for many purposes in our studies of electron motion, we do not have any 


great amount of margin, and that it will be necessary to use relativist 
dynamics whenever we wish to obtain really good results concerning 1! 
motion of even moderately high speed electrons. 

This article is purely expository. Its purpose is to set forth the funda- 
mental equations and theorems of relativistic particle dynamics in a cleat 


and concise form, unencumbered with any material relating to the theory of 
relativity proper. Almost all of the material is to be regarded as ety 
known, but apparently it is only to be found in an inconvenient and scattered 
form. The incomplete bibliography at the end of the paper gives references 
to some of the more accessible sources of this and other related material. 


IT. THe ELEMENTARY DIFFERENTIAL EQUATIONS OF MOTION 


Our discussion might be begun in any one of a number of ways, and no 
doubt the different approaches would appeal unequally to different readers. 
Considering the nature and purposes of this article, the author has deemed 
it best to begin by writing down at once the differential equations of motion 
of a particle (according to relativistic dynamics) in their most elementary 
form. Then, for the purposes of this discussion, these equations will have 
the status of a fundamental assumption. It need hardly be said that the 
equations are not written down arbitrarily. On the contrary, they represent 
the consensus of modern opinion as to the laws under which particles reall) 
do move.’ The grounds, experimental! and theoretical, for this opinion are 
set forth in various of the works cited in the bibliography. 

For the time being, until the contrary is stated in the final section, we 
employ a fixed rectangular coordinate system. Instead of denoting the 


coordinates of the particle by x, y, and z, as we have done provisionally 
the Introduction, we shall denote them by x, 22, and x3. Then 4, 4, and 
#3 denote the components of the velocity of the particle. The components 
of the force acting on the particle will be denoted by X,, Xs, and X3. For 
the time being we need only note that the force may depend upon the 


coordinates, the velocity, and the time; later on we shall introduce some 
more explicit assumptions about the force. The symbol ¢ will be used to 
denote the speed of light in vacuo. 

3 The validity of these laws is not unrestricted. It is limited on the one hand by the 
quantum phenomena which become appreciable on the atomic scale, and on the other hand 


by certain phenomena revealed by the general theory of relativity which becom: 
appreciable on the cosmic scale. 
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EQUATIONS OF ELECTRON MOTION 


We assume that the particle moves, under the influence of the force (.X, 


Y,, X3), so that its coordinates satisfy the system of differential equations 
d Mo Xn 


= X,, (wn = 1, 2, 3), (1) 
dty/1 — (v?/c?) 


where mo is a positive constant characteristic of the particle, and v* is an 


f abbreviation for the expression 7," + Z2” + 2@3;°.* The positive value of the 


square root is the significant one; and wherever square roots appear in the 
subsequent work it will be understood, unless the contrary is stated, that 
the positive values are intended. 

A few remarks may help bring out the significance of the foregoing assump- 
tion and its relations to the corresponding fundamental assumption of 
Newtonian dynamics. 

We call the constant mo the rest-mass of the particle, and we assume (in 
accordance with the experimental evidence) that mo is identical with the 
mass of the particle which is used in Newtonian dynamics. In relativistic 
dynamics the quantity m defined by the equation 


mo 
m= 


V1 — (e%/c?) 
is called the mass of the particle. We note that as v/c approaches zero the 
mass approaches the rest-mass (whence the appropriateness of the latter 
term), and that as v/c approaches unity the mass increases without limit. 

Consider the vector having the components 1, po, ps defined by the 
formulae 

a | Mo Xn 
V1 — (v*/c?) 
We call this vector the momentum of the particle. The momentum is equal 
to the velocity of the particle multiplied by the mass. 

Now equations (1) assert that the time-rate of change of the momentum 
of the particle is equal to the applied force. 

We have already observed that as v/c approaches zero the relativistic mass 
of a particle approaches the Newtonian mass. We now note that as 0/c 
approaches zero the components of the relativistic momentum approach 
the values 


Pu = mo dn (2') 
* We might merely say that v is the speed of the particle. However, for our immediate 
g \ 5 I I 


purposes, it is important not to lose sight of the fact that v is a certain particular function 
of the components of velocity. 
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which are precisely the components of the momentum according to 
Newtonian theory. 

Finally, as v/c approaches zero, the differential equations of motion | 
approach the forms* 


zi (mo En) = Xn; 1’ 
which are the Newtonian differential equations of motion. 

Thus we see that Newtonian dynamics is in effect a simplified approximat: 
form of relativistic dynamics which is valid when the speed of the part 
under consideration is sufficiently small compared with the speed of lig 

Let us carry out the indicated differentiations in equations (1), and the: 
solve the resulting equations for the quantities mod',, mod's, mod's. The work 
is straightforward, and need not be given here. We obtain the following set 
of formulae: 


xX, ii iec iya3¢° 
mit = (1-0)? | X, 1 — (a +47 dodo” 
X; tetsc 1— (ar +a )c° 
1— (a2 +4;)c° X; tidsc 
Mio = (1 - vc > ” 2, Toc . X» Xo X3C ’ " 3 
tidgc Xs 1— (a + &)c" 
1 — (¥-" + a3 )c , tidec’ Xx 
mis = (1 — vC = ss I) Lec . 1 — (x; -- a3 )c , Xo 
t,43¢°° det3¢° a 


These equations are, of course, the differential equations of motion (1 
written in a new, but equivalent, form. 
If, at some particular instant, the particle is moving parallel to the x-axis. 
so that #2 = £3; = 0, the equations (3) reduce at that instant to the forms: 
Mo X mo X2 . Mo X3 


a X - = X>, ——~. = X3. 
(1 — ve-?)3/2 1 (1 — vc—2)12 (1 — v%c-2)12 


These equations show that a particle of rest-mass mo, moving with speed 2, 


responds to a force parallel to the velocity as would a Newtonian particle 
of mass 
Mo 
m, = 
4 (1 — v%c-?)8/2’ 
4 If this conclusion is not entirely evident, the reader is referred to equations (3), fr: 


which the conclusion follows at once. 
5 T.e. an ideal particle which obeys the laws of Newtonian dynar irs. 
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d that the particle responds to a force perpendicular to the velocity as 
ould a Newtonian particle of mass 


m = 


For this reason, it was usual in the early work on relativistic dynamics to 
ascribe two masses to a particle: the ongitudinal mass m,, and the transverse 
mass m,. However, in general this procedure leads only to inconveniences, 
and it has been almost entirely abandoned. 

This concludes our discussion of the elementary differential equations of 
motion. Without any further general theory of relativistic dynamics it is 
possible to solve many interesting and important problems. For instance, 
it can be shown easily that the trajectory of a particle subjected to a force 
which is constant in magnitude and direction is a catenary (rather than a 
parabola, which is the curve predicted by Newtonian dynamics). In the 
following sections we shall discuss some of the less elementary parts of the 
subject. 


II]. THe LAGRANGIAN EQUATIONS 


In the foregoing the components of the applied force have been any func- 
tions of the coordinates, the components of the velocity, and the time. 
However, in problems concerning the motion of electrons, and for that 
matter in many other physical problems also, we are usually concerned with 
forces of a somewhat special kind. Throughout the remainder of the article 
we shall assume that the force belongs to this special class. 

We consider four given functions of the coordinates and time, namely 


V (x1, Xe, X3, t), A,(X1, Xo, ¥3, t), (72 = 2 ae 


and we assume that the components of the force are given by the formulae 


. OV OA, 5 OA» OA, ; OA, OA; 
Xy = —o — S84 gp) 2 — 1 a | SS 
OX) ot OX) OX2 OX: OX, | 
. OV OAs» _ | OAg Ole OA» aA, | 
X, = -— — — + hl ; —s ma Die = ’ (4) 
OX ot OX» OX | OX, OX2 | 
* AV ORs. OA, OAs OAs OAs | 
X3; = —- rel ii et) ae 7 & = 2583 = Ss 
OX2 al OX OX, OX» OX3 | 





Let us suppose, for purposes of illustration, that we are considering the 
motion of an electron. Then the physical interpretation of our assumption 


®*L. A. MacColl, American Mathematical Monthly, Vol. 45 (1938), pp. 669-676 
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concerning the force is the following. V(x, x2, x3, ¢) is the potential energy 
of the electron in an electromagnetic field; that is 


V(x, %2, x3, 2) = — ep(x1, Xe, X3, 2), 


where ¢ is the absolute value of the electronic charge, and g(x, ‘2, x3, ¢) is 
the scalar potential of the field. The functions A,(%, x2, x3, ¢) 2e related 
to the components a,,(%1, %2, %3, ¢) of the vector potential of the field by the 
equations 


A,(X1, Xe, X3, t) = €4,(X4, X2, X3, t). 


The terms —0A,/0t are —€ times the contributions of the vector potential 
to the components of the electric force. The quantity 0A3/0x_. — 0A2/d:; 
is —¢B,, where B, is the x,-component of the magnetic induction; and 
similarly for the quantities 0A;/0x3 — 0A3/dx, and 0A2/dx, — 0Aj/0.%.* 
In other cases also, equations (4), which may degenerate considerably, can 
be interpreted without difficulty. 
Now we define a function L(a, x2 


“) 


3, 4, Zo, £3, t) of the coordinates, the 
components of the velocity, and the time, as follows: 


L = — moc?(1 — v%c~?)'”® — V + aA, + bode + 29 A3. 5 


We call this the Lagrangian function. 
We write the equations 


d OL OL 


a 0 * = 1, 2, 3), (€ 
dt 02, OX n . ( » 4%, 9) 


carry out the indicated differentiations, and readily verify that the resulting 
equations are identical with those obtained by substituting the expressions 
(4) in equations (1). Hence, equations (6) are merely a form of the differ- 
ential equations of motion. We call equations (6) the Lagrangian equations. 
The chief importance of these equations is due to the ease with which they 
enable us to use coordinate systems which are not rectangular. This will 
be discussed in the final section. 

In the Newtonian case, i.e. the case in which the speed of the particle is 
small compared with the speed of light, the Lagrangian function reduces 
approximately to the form 


L= —mc ss ae (ay + a, + as) — V + & Ai + % Az + 43 43. (5’) 


* These relations between the A’s and the components of the vector potential, and 
between the partial derivatives of the A’s and the components of the magnetic induction, 
are based upon the use of the M.K.S. system of units. If we measure the electromagnetic 
quantities in other units, certain constant proportionality factors may appear in the 
relations. 
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If we employ the function (5’) in equations (6), we do indeed get the New- 


tonian differential equations. Since the constant term —mgc’ is of no 
effect in the formation of the differe:tial equations of motion, it is ordinarily 
omitted n writing the Newtonian form of the Lagrangian function. 

IV. HAMILTON’S CANONICAL EQUATIONS 


Let us write 


ba + An = Fn. 7) 





Solvii_, equations (2) for 41, %, #3, we get the result 


Cha [mo c+ pr 2 pe = ps | - 


Il 


In 


(ta — Aa)[me c+ (41 — Ar) + (42 — As) + (23 — As) T”. 


Also, it is readily seen that the differential equations (1) can be written, 
with the aid of equations (7) and (8). in the form 


OV o OA, . OAo ‘ OA; 
m=—- th. a2 ; t+ ds 5 
OXy, Ox n OX», Ox n 
(9) 
OV 


0 22 2 2 ‘ 2)1/2 
=— — c-—[moc + (m1 — Ar) + (m2 — Ao) + (93 — A3))] 
OX», 


Now let us define a function H(%), x2, x3, 71, 72, 73, t) as follows: 
H = clmoc’ +- (m - Ay ot (1 i Ao) ok (73 = A;)']' , + V. (10) 


Then equations (8) take the forms 


7 OH 
in = an,’ (11) 
and equations (9) take the forms 
‘ IH 
a = aot e (12) 
OX», 


The function H is called the Hamiltonian function. The six equations 
(11) and (12), which are equivalent to the three equations (1), are called 
Hamilton’s canonical equations of motion. These equations are of great 
importance in all of the deeper theoretical work in dynamics. 

An easy calculation shows that we have the identity 


H + L = Ww, 2, + 1222 + 13 X3. 
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In the Newtonian case the Hamiltonian function given by (10) reduces 

approximately to the form 
2 1 2 2 \2 , 

H = mc + = [(@1 — Ar) + (42 — Aa) + (43s — As) | + V. (10 
The equations (11) and (12), with H given by (10’), are equivalent to the 
Newtonian differential equations of motion (1’). Here again the constant 
term moc? is of no effect, and it is ordinarily omitted in writing the New- 
tonian form of the function H. The Newtonian forms of the functions // 
and L satisfy the identity (13), whether or not the constant terms my 
— myc are included. 


V. Static FIELDS OF Force: THE ENERGY INTEGRAL; NATURAL FAMILIES 
OF TRAJECTORIES 


By equations (11) and (12), we have the relation 


“ 3 ” 4 
dH _ 0H, ys K +. 4 oH | 


dt a «(A= Loz, Om, 
_ oH [aH dH dH dH] aH ‘s 
ot a=! 1 OL, Ofna OT, OXy Ot” 


In particular, if no one of the functions V, Ai, A2, As involves the time 
explicity, we have dH/dt = 0, so that the value of H remains constant 
during the motion of the particle. That is, under the condition stated we 
have 


me [1 — vc |” + Vix, x2, x3) = constant. (15 


In the Newtonian case equation (15) reduces approximately to the forn 
2 mo 2 , 
mc +—v +) (a1, %2, 43) = constant, 
which is equivalent to the equation 


mo 2 , sai 
7+ V (x1, %2, x3) = constant. (15 
It is well known that this equation is a consequence of the Newtonia! 
differential equations of motion. 
. The left-hand member of equation (15’) is the energy of the particle 
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Newtonian dynamics, the first and second terms being the kinetic energy 


and the potential energy, respectively. The equation itself is called the 
energy integral.’ Similarly, we call (15) the energy integral in relativisti 
dynamics, and we call the expression 

mec({(l—wve |] “+ V 
the relativistic energy. This energy is the sum of three parts: the proper 
energy moc”, the relativistic kinetic energy 


mo [1 —ve |)” — mo 


and the potential energy V. 

The totality of possible trajectories of a particle in a static field of force 
forms a five-parameter family. We now see that if the field of force is 
static and of the kind we are considering now, the five-parameter family of 
curves consists of « ! four-parameter subfamilies, each of which corresponds 
to a different value of the energy of the particle. Each of these four-param- 
eter subfamilies is called a natural family of trajectories. We proceed 
to derive the differential equations defining a natural family. 

If the constant in the right-hand member of equation (15) is denoted by 
the symbol £, we have the relation 


x,{1 + Xx 


42 
“» 


+ x3]? = ct — mct(E — V)°}", 16) 
where 

x2 = dx2/dx, x3 = dxs/dx. 
Hence, 


dt=c(itxe +23] 71 — mic (E-— V) 7)” dx. 


From this, and the two equations 


d Mo Xo OV a OA3z OA» , OAs OA, 
“aoe OC ieee ales wate ar se hae eae | fay 4 

dt (1 — we-?)'? OX» OX2 OX3 Ox) OX 

d Mo X3 OV : OA, OA3 : OAs; OAg 
; Sees, ae: ee eat UY = = — ae a es 

dt (1 — 9c *) 1/2 OX: OX3 Ox, OX» OX3 


it follows that we have the following system of differential equations defin- 
ing the natural family of trajectories corresponding to the total energy E: 


7In the theory of differential equations, an equation relating the unknowns involved 
na system of differential equations, their derivatives of orders less than the highest orders 
appearing in the system, the independent variable, and one or more arbitrary constants, 
s called an integral of the system of differential equations. 
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ee ee ee 
te +e TS («:| cf Sida 
dx 1 + X2 + X3 
: ; y 
= 5 (E — Vy = mict}"* 
Ng #2)—1/2 uy dA» VAg 
T cll + te + &s | = («i oA - es le one ~ Is = OA, ), 
OX OX3 Ox, OX» 
2 r2 ea ’ Pom 2 2 2s 
(1 + Xe + X3 ia wes (2: k = | ) 
dx lta +45 


0 ‘ ‘ m 
= —[(E — V) — moc'|'” 
on, | ) oc | 


2 ey i, wy ' A; VA» 
+cit+ Page oe ‘fie (| on : J — X2 k ae 8 )). 
Ox Ox, OX» OX3 


The equations which correspond to (17) in the Newtonian case are most 


li 


readily obtained by going back to the Newtonian differential equations of 


motion and employing the integral 


mov?/2 + V = I 


An easy calculation, which is entirely parallel to the foregoing, gives us the 


following system of equations: 


tas + To (|. : ani \=% (E - vy” 
dx, 1+ i ~f- x3 Ax. 
mo) , ig re fe: 1/2 Bg OA; = OA» ee OA, ) 
+ [2mo(1 + 12 + 13) € [fe 24s] a OX» : 17’ 
(1 / 
(1 + x2. +e7° CF x »| )- © (EB - v)"" 
= 1+ xo z xs OX 


12 eet OA 0A 110A OA 
! J 1/2 eos i. a 2 
+ [2mo(l + x2 + a )] (| 0X3 =] = E “asl. 


On comparing the systems of equations (17) and (17’), we get the follow- 


ing useful theorem. 
If the constants E, E*, mo, mo*, k, and the functions (of x, X2, X3,) V, A 
As, V*, Ay*, Ao*, A3* are such that we have identically 
’ ’ - 


(E—V) — mci = R(E* — V*), 


04; 942_ sik 04; dA? 
OX2 0x3  -c(2mo*)"/2 | Axe Ox3 |’ 
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a o e — e 

dA, OA, k 0A> = AAT 
Ox, OX» c(2mo*)"/? | Oxy OX | . 
« - « + ¢ oa 
OA, OA3 7” k OA, OA; 
Ox Ox, c(2mo*)'/? | Axes Ox, |’ 


the natural family of trajectories of a relativistic particle’ (of rest-mass m 
moving with relativistic total energy E in the field of force derived from the 
functions V, Ay, Aa, Ag is identical with the natural family of trajectories of a 
Newtonian particle (of mass mo*) moving with Newtonian energy E* in the 
field of force derived from the functions V*, A,*, A.*, A;*. 


In particular, the conditions of the theorem are satisfied if 


1/2 . ‘ o 2 4 
k = c(2mo) ”, E* = c-2(2mo)—(E2 — moc), mo* = Mo. 


Hence, we have the corollary: 

In the case of an electrified particle moving in any static magnetic field the 
natural family of trajectories corresponding to any value of the energy given by 
relativistic (Newtonian) dynamics is identical with the natural family of 
trajectories corresponding to a certain other value of the energy given by New- 
tonian (relativistic) dynamics. 

The equation 

E* = ¢*(2m) (E’ — moc’ 
establishes a one-to-one correspondence between the physically significant 
E = moc? and E* = 0) values of the relativistic energy E and the New- 
tonian energy E*. From this fact and the preceding corollary we get the 
following further result: 

In the case of an electrified particle moving in any static magnetic field the 
total five-parameter family of trajectories given by relativistic dynamics is 
identical with that given by Newtonian dynamics. 

Of course, these peculiar properties of motion of an electrified particle 
moving in a static magnetic field are explained physically by the fact that 
the magnetic forces do no work, so that the speed of the particle, and 
consequently also its mass, remain constant during the motion. 


VI. Some FORMULAE FROM THE CALCULUS OF VARIATIONS 


This section is devoted to the derivation of some formulae from the 
Calculus of Variations which will be needed in the further discussion of the 


*T.e. a particle obeying the laws of relativistic dynamics 
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dynamics of a particle. All constants, variables, and functions considered 
here are understood to be real.” 

Let F(t, x, y, 2, p, g, 7) be a function of the seven arguments indicated. 
which, together with all of its partial derivatives of the first three orders, ; 
continuous in a region R defined as follows: 


a4<t< a, 
by << t < de, 
Ri a<y<@, 
a<z< dh, 
p, 7, and r unrestricted, 


the a’s, b’s, c’s, and d’s, being constants. 

Let x(t), y(t), s(t), o(t), W(t), and w(t) be continuous functions with con- 
tinuous first derivatives, and let €, 7, and 6 be parameters, independent of 
such that we have the relations 


by < x(t) + elt) < be, 
a < y(t) + mt) < ©, (aq) < t < dy). 
d, < s(t) + O@w(t) < do, 
Let 7, and 7, be constants, and let 4; and ¢, be parameters, such that 
a4<7Ty,+h < To + be < &. 
We now consider the integral 


Te, Un 6, hh, to) 


Totte 
= i F(t, x + eg, y + nv, 3 + Ow, x’ + eg’, y’ + YW’, 2 + Ow’) di 
T 


ratty e 


It can be shown without difficulty that the integral exists and is a differen- 
tiable function of €, 7, 6, ti, fz. We are interested in formulae giving the 
values of 07/de, 01 /dn, 01/00, AT /dt;, OT /dt, at the point € = n = 6 = ft, = 
=> Q. 

® Since this section is purely mathematical, the constants, variables, and functions 


not necessarily have any special physical significance. 
1° We treat the case of a function of seven arguments in order to fix the ideas, and 


because this is a case we shall meet in Section VII. However, the discussion applies 
essentially to other cases as well. In particular, in Section VII we shall also deal with a 


case in which F has only five arguments, z and r being absent. 
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By a well known theorem concerning the differentiation of definite 
integrals with respect to parameters,'! we have 


‘ Tette 4 
ol i d , 0 , ' , 

= ¢ ~+¢ - ~ | F(t,x + eg, ---, 2’ + Ow’) di, 
de Ty+t} [ O(x + €¢) ” O(x’ + ey) . ™ 


: = —F[T, +4,2(71 +h) +e¢(7T1 +h), +++, 2'(Ti +h) + Ow'(T + 4)], 
hy 
t = F(T; + bh, x(T2 + bh) + e9(T2 + bh), +++ (2 + te) + Ow’ (Te + be). 
ty 


The formulae for 0J/dn and 07/06 are similar to that for 07/0, and need 
not be written down. 

In particular, if [87/de]o, etc. denote the values of the derivatives at the 
pointe= n= 0=%4 = kh = 0, we have 





ar | a 0 0 
= A ,/ Fit, x, ---, 3') dt. 
rt 4 E ax +¢ | (é, 2, It 
al) "2 a) ) 
= | lv —+W - ] F(t, x, «+ 2!) dt, 
| On_jo qT; oy oy 
ar] f° te) >, @ ; mans 
| 06 _jo a3 T; E Oz =. 3] pets 9 SR - 
_ 0 = —F\T,, x(T)), --- ,2 (Ty), 
L778 
Pas” 
oe F([T2, «(T2), «++ , 2’(T2)). 
Lr 








The first three of equations (18) can be transformed to advantage, as 
follows. Integrating by parts, we obtain the formula 


T, F) T2 
a! ; - e oe! = ate eee mad 
| ¥ ax’ F(t, x, a )dt & Ox’ PF (t, Ys i 


- I". i= Ps a, 
and similar formulae for the integrals 


i. VP: ak 


TT? F) 
i a! © F(t x, --- ,#) dt. 
7 Oz 


dl 


nn 


and 


'' The theorem is given, often in the form of two separate theorems, in most works on 
\dvanced Calculus and the Theory of Functions of Real Variables. See the bibliography. 
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It follows, therefore, that we have 


0 ‘y T2 
|] we E he , F(t, x, °**, ‘| 
de jo Ox . 
T? 
/ 0 d a ; 
-+ a(t — wecasn F(t, x gt ia 
I o ) K dt 2] x, ) dt, 
7 ee Te 
| = lv o. F(t, x, «°° : “| 
On _jo oy rn 
“wll : : F(t, d 
Ty dy dt dy’ a) »2) at, (19 
a) ) Ts 
dé jo Oz i 
T2 4 * 
_ ce) doa : : 
: I wo) E dt 3 | P(t, x, +++ 8’) dt. 


An important special case is that in which 4, and f are zero (so that th: 
limits of integration are fixed), and 
e711) = (72) = (Ti) = (72) = w(T)) = w(T2) = 0. 


In this case we have in general 


/ OF d OF 
I(e, n, 6,0, 0) — 1(0, 0,0, 0,0) = fe oF _d . 
bicchaaiaate, Ty w | dt ax’ | ‘ 


OF d aA F 
lis ‘s vi {or - 4 dt dy =| + of “7 | - dt 02 ‘| 


+ o(e, n, 4), 


where o(e€, n, 8) denotes a term, the exact form of which is unimportant 
which is such that the expression 


o(e, n, 8) 
le] +] n| + |4| 


approaches the limit zero as €, n, and 6 tend simultaneously toward sero. 


In particular, if the functions x(t), y(t), s(¢) satisfy the system of differen- 


tial equations 


d0F_ OF _, d@9F OF _, daF_ oF 


a en ee, = 0, me (0 
dt dx’ Ox dt dy’ oy dt dz’ 02 ; 


we have (for all choices of the functions ¢, ¥, w subject to the conditions 
stated) 


I(e, n, 8,0, 0) — I(0, 0, 0, 0, 0) = ofe, n, 8). (21 
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Also, it can be shown without difficulty that in order that we have (21), for 


all such choices of ¢, y, w, it is necessary that x, y, and s, satisfy the equations 
0).* 


The last result can be stated in the following summary, and not quite 


F explicit, form: If, and only if, the functions x(¢), y(t), s(¢) satisfy equations 


20), the integral 


T? 
| F(t, x, y, 2, x’, y’, 2’) dt (22) 
Ty 
is stationary with respect to infinitesimal variations of the functions x(¢), 
vif). s(t) which leave the terminal values unaltered. 

The problem of finding functions which render the values of definite 
integrals stationary is the chief subject of the Calculus of Variations. 

The equations (20) are called the Eulerian equations of the Calculus of 
Variations problem of making the value of the integral (22) stationary, or, 
as we usually say, of maximizing or minimizing the integral. 


VII. HAMILTON’S PRINCIPLE AND THE PRINCIPLE OF LEAST ACTION 


We immediately recognize equations (6) as the Eulerian equations of a 
problem in the Calculus of Variations. Thus we have the following principle 
Hamilton’s principle): 

The particle moves, under forces of the type (4), so that the value of the integral 


te 
| L dt, 


with ty and ty held fixed, is stationary with respect to infinitesimal variations of 
the functions x,(t) which leave the initial and final points unaltered. 

The precise meaning of this is determined by the discussion given in 
Section VI. 

Hamilton’s principle leads to the relativistic or Newtonian differential 
equations of motion, according as we use in it the function Z given by (5 
or by (5’). 

A little inspection suffices to show that the system of equations (17) is 
also the system of Eulerian equations of a problem in the Calculus of 

ddF OF 


* In brief, suppose that —. — —— were not zero for some value of ¢. Then if we should 
dt ax Ox 


choose a function g(t) which was (say) positive in the neighborhood of that value, and zero 


elsewhere, the integral 
i ) dF d@F 
t _ toll 
a . Ox dt 0x’ 
1 


would have a value other than zero. We shall not give the actual proof here; it is to be 
found in the works on the Calculus of Variations cited in the bibliography. 
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Variations. Thus we get the so-called principle of least action, which 
stated as follows: 

The particle moves, in a static field of force of the type (4), and wit) 
prescribed total energy E, in such a curve that the value of the integral 


where 
functic 
It fc 
F of the 
(2 do 
a 21/27 r\2 —2 2 2)1/2 4 ' 
| (1 +a. +43] ((E—V)'c° — moc) + Ai + Aoxe + A325) dy L 
(41), 
with the limits of integration held fixed, is stationary with res pect to infinitesim 
variations of the trajectory which leave the end points unaltered. 
We have a precisely similar principle in Newtonian dynamics, but her 
the integral in question is 


(#1) 2 42 > 2 ’ , 
| : ({1 + Xe + X3 |! “(2mo(E stad V)|' . ca Ay fe Ao Xe a A3%3) dx j 
( 


71), 


The last two integrals can be written more symmetrically, but not quite 


so explicitly, as follows: 


Pe 
. »\2 —2 2 21/2 1: 12» d 
[ (1 - VV — mic}! + AB + 4488 + 4,8) as 


i 


Po 
. » 2 d x d. Xe dx: 
| [2mo(E — | )}" + A; — + Ao— + Az — ds, 
Py ds ds ds : h: 
so thi 

where P; and P, denote the end points of the trajectory, and ds? = d i+ point 
dx? + dx3. are t 
form 

VIII. THe HAMILTON-JACoBI THEORY 


Let us write 


W = i ~ Lai (t), xa(t), x9(t), xi(2), xa(t), x4(0), tl dt. 


t 


We have already studied the variation of W when ¢; and f& are held fixed 

and the functions x,,(¢) are varied in such a way that the terminal values are 

unaltered; and we have shown that under these circumstances the variatio! 

of W vanishes, to the first order of small quantities, in the natural motion." whic 

In the following we shall study the variation of W under some other & AW 

conditions. 
Specifically, we shall study the quantity AW defined by equation (23 

and the equation 


tatAte 


W+ AW = Liai(t) + &:(t), «++ , x(t) + &3(2), €] dt, 


ty+St, 


2 T.e. a motion satisfying equations (1). 
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where the functions x,(¢) represent a natural motion, the &,(¢) are small 


functions, and Aé, and Af are small parameters. 
It follows from the results of Section VI that we have (to within terms 


of the second order in small quantities)'* 
Abs Lx; (te), ek x3(to), to] 


— At, L|x(t), «+: , x3(t), ty] 
3 : 
ii E(t) 
t=t 


3 
OL 
+ 2 fell lo ; n=] Ox’ 


n=l 
Ab Lixi(t2), «°° , x3(t2) to] 
— At L[x(t), «++ , x3(h), tl 


3 
+ bs [an (te)En (te) ee mn (tyEn(tr)| 
n=l 
Let us write 
(Arno = Xnlbe + Ats) + E, (le + Ate) — X,( le) nits) + xa (te) Ab, 
(Ax,)1 = Xnlh a At) 4 En(h + Ad, ) — 2X, (h) 
so that (Aa,)o, (Axe), (Ax;)2 are the coordinate differences of the terminal 


, 
En(t1) + Xalh )At, 


points of the varied and unvaried curves, and similarly (Ax,);, (Ax2):, (Ax3); 
Then we have the 


ire the coordinate differences of the initial points. 


formula 
: , 
AW = (tints), J > ra(ts)2n(t)) Ale 

n=) 


3 
-J- > ra(ti)3s(h)) At, 


_ (din(n), 
n=l 


3 
+ DS [an(te)(Axn)o — ta(tr)(Axn):] 
n=l 


which, by equation (13), can be written in the form 


—H|x;(t2), so -] Ab + H (x(t), ig -] At 
I. (24) 


AW = 
3 
+ > [aa(te)(Axn)2 — n(tr)(Axn):] 


n=l 


Now, the integration in (23) being taken over a natural motion of the 
particle, the value of W depends upon the initial instant, the initial coordi 


'’ This is also the sense in which the following equations are to be understood. 
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nates, the initial components of velocity, and the final instant. It 
necessary now to consider W as depending upon the following equivalent se 
of eight variables: the initial and final instants 4; and f, the coordinates 
(X11, X21, Xs) of the initial point, and the coordinates (12, 422, X32) of the ti 
point. Regarding W in this manner, we at once obtain the fo 
relations from equation (24) 


ow Ow 

. = — Ho, ; = Tn2, 2) 
Oly OX n2 

ow Ow 

- _ Hh, ™ = —Tniy, yf 
Ol, OXn1 


where H, denotes H{x;(t2), «++ , m(te), «++ , | and H, denotes H[x,(¢,), 
Ml), °°*, 4]. 
Let us now consider the partial differential equation 

ow eee Sponge rr : 
ot + I(x, » 2, X, ou Ox, ov / OX2 , ov /OX3 ,4= Q). 2 

€ 
The preceding work shows that the function W we have been considering 
(with an, %21, Xs, 4 regarded as parameters, and with the symbols x», 


X32, fg replaced by x1, %2, x3, ¢ respectively) is a particular solution of this 


equation. We shall show that the complete solution of this equati 
possesses remarkable properties in connection with dynamical problems 
The complete solution of equation (27) is a function of x, x2, xs, f, al 
of four arbitrary constants, of which one is merely additive, and can |x 
neglected for our purposes. Let the solution be written 
W = Wa, Xo, X3, l, a1, Qe, a3), 
where the a’s are the three essential arbitrary constants. 
We write the equations 
aw ) 
= Bn, (<5 
Jan, 


‘ 


where the @’s are further arbitrary constants. These equations implicit! 


determine the x’s as functions of ¢ and the six arbitrary constants ay, - - 
We also write the equations 
OW _ 
OXy 


oT 
ae (29 


These equations determine three functions 7, of the «’s, the a’s, and / 


In virtue of equations (28), the z’s are ultimately functions of ¢, the a’s, a! 


the B's. 
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[here is no reason to foresee a priori that the functions x(t, a1, --- , 83), 
. , wa(t, a1, *** , 83) determined in this way, by means of the complete 
solution of equation (27), satisfy the differential equations of motion (11 
and (12). Nevertheless, they actually do satisfy those equations, as we 
proceed to show. 
By equations (28), we have the relations 
dB, oW .. ow : 


0 = = —— ae ae W 
dt Oa, ot mal On OXms 


- 14 
On the other hand, by (27) and (29), we have 


0 OW 
GQ = . ' - + H(x,, X2, X3, Wi, We, 73, | 


Oan ol 

a2 . 3 o ” «2 , 3 e alrrr 

ou OH Om ow , > oH av 1 
= - = + . (3 

0a, Ot wnt Olm Olle da, Ol ment OR mn Olly Ohm 


The determinant 


aw aw 
day Ox 1 day OX 

ow ow 
Oa OX, Jag OX3 


isnot zero. For if it were, we would have a relation of the form 


~ 
i) 


PAW /dx,, OW /dx2, OW /Ax3, x1, X2, Xs, t] = O, 


independent of the a’s. Now equation (32) is obviously distinct from (27), 


since it does not involve dW /dt. Hence, the vanishing of the determinant 
would imply that the function W (a, x2, as, t, a1, a@2, a3) satisfies two distinct 
partial differential equations of the first order. This, however, is impossible 
when W is the complete solution of (27); for an essential part of the concept 
of the complete solution of a differential equation is that the elimination of 
the arbitrary constants, from the solution and the equations obtained by 
differentiation, shall result in the given differential equation and no other. 
It follows, therefore, from (30) and (31) that 


OH 
‘ a 
We also have, by (29), 
; aw ow OW , ow -. #W OH a 
mm Fe ~~ d = -; Py - - s (35 
¥ Ox, Ot iz = nt, Ox, Ot m=1 OXm OXn OF m 


4 Since the function W (x1, Xo, X3, t, aj, a2, a3) satisfies equation 27) identically in the 


v's, t,and the a’s. This remark applies also in the case of equation (34) 
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On the other hand, we have 


o a » a2 , o 3 e Quy, 
leas & +1] _ ol , oH Me.» 0H oN : 


Ox, | Ot Ox, OF OX, ml O8m OX OL 


By (33) and (34), we have the second set of canonical equations 


i- OH 
Ox, 


™ = 
This completes the demonstration. 
If H does not involve the time explicitly, we can write 


W=S- Et, 3s 


where £ is an arbitrary parameter, and S is a solution of the differentia 
equation 


H[x1, Xo, X3, Os OX, Os OX, Os 0x3] = Ff. 36 


The complete solution of (36) contains three arbitrary constants (besides 
the parameter £), of which one is merely additive, and can be neglected 
It is easily seen that the solution of the canonical equations determined i: 
the way described above, using the function W given by (35), and treating / 
as one of the a’s, represents a motion of the particle with the total energy 

All of this theory holds both for the relativistic case and for the New 
tonian case, the only difference being in the forms of the differentia 
equations (27) and (36) in the two cases. 


LX. CURVILINEAR COORDINATES 


In all of the foregoing we have employed rectangular coordinates, because 
they afford the simplest and most direct expression of the basic physica 
facts. However, in the solution of particular problems it is often more 
convenient to use other systems of coordinates. For this reason, we shal 
now formulate the more important equations in terms of general curvilinear 
coordinates. In this work, as in all work with general coordinate systems 
we shall encounter concepts and relations which can be handled most 
perspicuously by means of the modern tensor calculus. Actually, the 
amount of tensor calculus we shall use is very slight, and no extended pre- 
liminary discussion is necessary in order to make the formulae intelligible 
It will suffice to give occasional explanations of the notation, and of som 
of the concepts, as we proceed. Further information is to be found in the 
works cited in the bibliography. 

First consider the Lagrangian equations, which, as we have seen, are 
merely the Eulerian equations which follow from Hamilton’s principle 
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Now Hamilton’s principle expresses a fact concerning the motion of a 
particle which is, by its very nature, independent of the choice of coordinates 
Hence, the Lagrangian equations (6) hold in any coordinate system. How- 
ever, the form of the function 1 depends upon the particular coordinate 
system, and we must discuss the change of the form of the function resulting 
from a transformation of the coordinate system. 

In accordance with the common practice in the tensor calculus, we shall 


° 1 2 a P 
now denote the coordinates by the symbols x, x°, x, instead of by the 
symbols 44, Xe, ¥3. 


In rectangular coordinates the differential distance ds between the points 


tx?) x3) and (x! + dx!, x? + dx?, x3 + dx) is given by the simple formula 
. 2 a? 2 
ds? = dx + dx? aa dx3 ‘ 


but this is highly special; in general coordinates we have 


3 3 
2 Fe 2 3 n - 
ds = pm be Smn\X > Xx ) dx” dx ; (37) 
m=1 n=l 
where the g’s are functions which depend upon the particular coordinate 
system under consideration. It is understood that gm, = Zam.  Hence- 
forth, we shall write (37) in the form 


ds? = g,,(x', x*, x*)dx"dx", 38) 
and we shall observe this general rule throughout: When the same literal 
index occurs twice in a term, once as a subscript and once as a superscript, 
that term is understood to be summed for the three values of the index. 

We now have the result 


o? = [ds/di|? = g..,(x', x?, x*)Z"2", 
and 


9 


moc?(1 — vc 2)1 3 as Mol “11 —C¢ ae ad 2 


The function V(x', x, x*, 4) isa scalar. That is to say, when the coordi- 
nate system is changed, the first three arguments of the function are replaced 
by their expressions in terms of the new coordinates, and so we obtain a 
function which is of a new analytical form, but which has the same value 
as the original function at each point of space. 

Now we consider the expression 


A,z! + Aox? + A 32°. 


In rectangular coordinates this is the scalar product of the vectors (A), 
A», As) and (z!, #*, 8). The expression retains its form and interpretation 
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under changes of the coordinate system, provi ided (as the notation implies 
(A, Ag, A3) is treated as a covariant vector.’ 

With these understandings as to the significance of the symbolism, we 
can now write down the following general expressions for the Lagrangia: 
function Z in the relativistic and Newtonian cases, respectively, 


L = —moc [1 — € gunk” z’|'” — V + Ant”, 


2 , Mo .m on , 
L=—mce + 5 fant £ — V+ A,” 


These hold for any coordinate system; and from the appropriate one of these, 
and the Lagrangian equations 


daL_ aL _, 
dt az” ax" ”’ 


we obtain the relativistic or Newtonian differential equations of motion in 
any coordinates. 

Now let us consider the Hamiltonian canonical equations. 

We have already agreed to consider (A;, Ao, As) as a covariant vector 
We now make the same convention in regard to (m1, #2, 73). Then it read- 
ily follows that the equations 

OL 
ain ™ (40) 

18 Suppose that with a point P (which may be either a special point or a typical point 

and with each coordinate system, we have associated an ordered triple of numbers. 


If the triples of number (a, @2, @3) and (a,’, @2’, a3’) associated, respectively, with an\ 
two coordinate systems (x, x?, x%) and x!’, 2?’, x3’) satisfy the relations 


Ox” 
Int = ax™’ Any 


the numbers (a), @2, a3) are said to be the components of a covariant vector in the coordinate 
system (x!, x?, x3), (It is understood, of course, that the partial derivatives are evaluated 
at the point P.) 

On the other hand, if the triples of numbers (a', a*, a*) and (a!’, a?’, a3’) associated with 
the typical coordinate systems (x!, x?, x) and (x!’ , x3’) satisfy the relations 


Ox™ 
am’ — a", 


the numbers (a!, a, a’) are said to be the components of a contravariant vector in the 
coordinate system (x!, x?, x). 

These concepts agree only in part with the ones used in the elementary theory of 
vectors. From our present standpoint, the only vectors used in the elementary theories 
are those which are defined with reference to rectangular coordinate systems. When other 
coordinate systems are used (e.g. cylindrical coordinates), the vectors, defined in terms of 
rectangular coordinates, are merely resolved along the tangents to the coordinate curves. 
The components obtained in this way are not the same as the components considered ir 
the tensor calculus, which we are using here. 
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are tensor equations; and since they hold when the coordinates are rec- 


tangular, they hold for all coordinate systems.” 
We let g™" denote g~! times the cofactor of the element g,,, in the deter- 


Sr 


minant 


Now we write 
H = cmc + g""(1m 
for the relativistic case, and 
H = moc? + (2mo)—!g"" (atm 
for the Newtonian case. We see that these expressions specialize into the 
ones given earlier for the Hamiltonian function when the coordinates are 


rectangular. 
H, L, and r,2" are all scalars. Consequently, the equation 


H+L= r,2" 42) 


is a tensor relation; and since it holds when the coordinates are rectangular, 
it holds for all coordinate systems. 
The Lagrangian equations can be written in the form 


= OL (43) 


Ox" ” 
Let us consider the variation of the function Z resulting from small 


variations of the x’s and 2’s._ By (40) and (43), we have the relation 


aL oe | iia 
gr ee gat a 


Ox" Ox” 
= 7, 6x" + 2,62" 
= §(x, 2") + (w,6x" — x” bn, ). 


It follows from (42) and (44) that the variation of H resulting from small 
variations of the x’s and the z’s is given by the formula 
5H = Zin, — 7,02". 


6 The argument is explained in detail in the works on the tensor calculus cited in the 
bibliography. 
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From this it follows that we have the Hamiltonian canonical equations 


ne oH F OH 


tm = — 
’ n ’ 
Onn Ox" 


in any coordinate system. 

We have already seen how to state Hamilton’s principle in terms of gener 
coordinates. 

In the relativistic case the principle of least action takes the form: 7); 
particle moves, in a static field of force of the type (4), and with the prescribed 
total energy E, in such a curve that the value of the integral 


me dx” dx" |! ‘ : -\2 2 2 21/2 dx™ ; 
zs (| ss dx' = \(E = J F< — Mo ( + An ir) a - 


with the limits of integration held fixed, is stationary with respect to infinitesimal 
variations of the trajectory which leave the end points unaltered. The corre- 
sponding form of the principle for the Newtonian case is obvious. 

We are now in a position to dispose very quickly of the problem of formu- 
lating the Hamilton-Jacobi theory in terms of general curvilinear 
coordinates. 

The general form of the Hamiltonian function being given by (41) (for 
the relativistic case) or (41’) (for the Newtonian case), we can at once write 
down the partial differential equation 


ow 
ot 


+ H(x', x’, x°, AW/dx', AW/dx’, AW/d2’, t) = 0. (45) 


Let 
ram Wat at ae i 
W = W(e', x3, x*, t, a’, a*, 2°) 


represent the complete solution of (45), without the irrelevant additive 
constant of integration. 
Our chief problem is that of proving that the functions a"(f, a!, a’, a’, 
81 , Be, 83), mrlt, a', a, a*, By , 82, B3) determined by the equations 
Ow _ ow 


7 Ba, — Ses 


da" Ox” 


where the ’s are further arbitrary constants, satisfy the canonical equations 


.» _ OH < OH 
a 


= Te ——. 
Or,’ Ox” 


Now, referring to the proof given in Section VIII for the special case of 
rectangular coordinates, we see at once that nothing in the proof depends 
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upon the special forms which the Hamiltonian function and equation (45 
assume in those coordinates. Hence the proof already given applies imme- 
diately to the present general case. 

Similar remarks apply also to the case in which H does not involve the 
time explicitly, and in which we write 


W=S-— Et, 


where S is the complete solution (without the additive arbitrary constant 
of the equation 


H(x', x7, x3, 0S Ox!, 0S/dx?, 0S/dx?) = E. 
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CHAPTER I 
Quartz Crystal Applications 
By W. P. MASON 


1.1. INTRODUCTION—PURPOSE OF SERIES 


HIS paper is the first one of a series of papers dealing with quartz 

crystals, their applications in oscillators, filters, and transducers, and 
the methods of producing them from the natural crystal. This series was 
prepared first to make available to the Western Electric Co. and other 
manufacturers of quartz crystals some of the specialized knowledge on 
these subjects that has been acquired at the Bell Telephone Laboratories 
Sufficient interest has been expressed in this series to make it desirable to 
publish them in serial form. 

This first paper in the series is a general introductory paper covering 
the application of crystals to oscillators, filters and transducers. An appen- 
dix is given which discusses the elastic and electric relations in crystals and 
gives recent measurements of the elastic constants, their temperature coeffi- 
cients, and the piezoelectric constants of quartz. This paper is followed 
by more detailed papers by Messrs. Bond, Willard, Sykes, McSkimin, and 
Fair which give consideration to quartz crystallography; determination of 
orientation by optical methods, X-ray methods, and etching methods; the im- 
perfections occurring in quartz crystals; modes of motion and their calcu- 
lation; the dimensioning of crystals to avoid undesirable resonances; and 
the use of crystals in oscillators. 


1.2 Earty History OF PIEZOELECTRICITY AND ITS APPLICATIONS 

The direct piezoelectric effect was discovered by the brothers Curie in 1880) 
They measured the effect first for a quartz crystal by putting a weight on the 
surface and measuring the charge appearing on the surface, the magnitude 
of which was proportional to the applied weight. A simple model for demon- 
strating this effect can be made by using a large piece of Rochelle salt cut 
with its length 45° from the Y and Z crystallographic axes and placing tinfoil 
electrodes normal to the X axis. If these electrodes are connected to a 
neon lamp, and the crystal is compressed by hitting it with a hammer, a 
charge is generated on the surface and a voltage applied to the lamp sufficient 
to break itdown. In fact as much as 2000 volts can be generated by striking 
the crystal hard. 
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The converse piezoelectric effect was predicted in 1881 by the French 
physicist Lippmann on the basis of the principle of conservation of elec- 
tricity. It was verified in the same year by the brothers Curie. In this 
effect a crystal is strained when a voltage is applied to it. The effect can 
be demonstrated by a model which consists of two thin pieces of Rochelle 
salt poled so that one expands when the voltage is applied and the other 
contracts. The result is—as in a bimetallic thermostat—the crystal bends. 
For crystals 10 mil inches thick and 4 inches long, a ninety-volt battery 
applied causes a displacement of a quarter of an inch or more of the end of 
the unit. Reversing the voltage reverses the direction of the displacement. 
The Curies constructed a bimorph unit of this type out of quartz and used it 
practically to measure voltage by measuring the displacement of the end 
of the crystal. By connecting the leads of an electrometer to the terminals, 
they could measure force applied by measuring the amount of charge gen- 
erated at the terminals. 

Outside of this use which was quite minor, the piezoelectric effect remained 
a scientific curiosity until the war of 1914-1918. It did inspire, however, 
considerable scientific speculation. Lord Kelvin in 1893 proposed a model 
for explaining the piezoelectricity of quartz and was able to calculate 
approximately the value of the piezoelectric constant. This model is 
discussed briefly in the next section. He also constructed and demonstrated 
a ‘piezoelectric pile’? made from small spheres of zinc and copper, to illus 
trate the effect. At about the same time (1890-1892) Voigt published a 
series of papers followed by a book “‘ Lehrbuch der Kristal! Physik” (1910 
in which the stresses, strains, fields and polarizations of piezoelectric crystals 
are related in mathematical form. These mathematical expressions (which 
are discussed further in the appendix) form a basis for the development of 
the properties of oriented crystals as discussed in section 1.5. 

During the war of 1914-1918, Professor Langevin in Paris was requested 
by the French Government to devise some way of detecting submarines 
by acoustic waves they produce in water. After trying several devices he 
finally found that piezoelectric quartz plates could be used for that purpose. 
His device, which is shown in Fig. 1.1, consisted essentially of a mosaic of 
quartz which has the property that when a voltage is applied the crystal 
will expand and send out a longitudinal wave. Similarly, if a wave strikes 


it, the wave will set the quartz in vibration and generate a voltage which 


can be detected by vacuum tube devices. Langevin did not get his device 
perfected till after the war so it was not used at that time to detect sub- 
marines. Similar devices have, however, been used in this war. Langevin’s 
original apparatus was used extensively as a sonic depth finder. In this 
use a pulse is generated which is recorded directly on a moving record and 
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is also sent out into the ocean. It strikes the bottom and is reflected back 
causing another mark to appear on the record. Knowing the difference 
time and the velocity of sound in sea water, the distance to the bottom ca; 
be measured. A typical record is shown in Fig. 1.2. The top record shows 
the contour of the sea bottom while the second record shows the reflections 
from a school of fish. 

At about the same time, Nicolson at Bell Telephone Laboratories was 


experimenting with Rochelle salt, another piezoelectric material having 





Fig. 1.1—Ultrasonic transmitting apparatus 


strated loud speakers, microphones, and phonograph pick-ups using Rochelle 
salt.'_ He was also the first one to control an oscillator by means of a crysta 


between the resonating coil in the output and the grid, since the electrode 


1“The Piezoelectric Effect in the Composite Rochelle Salt Crystal”—A. M. Nicolsot 
Proc. A. I. E. E. 1919, 38, 1315. 
2 See Patent 2,212,845 filed April 10, 1918; issued Aug. 27, 1940. 


much larger piezoelectric effect than quartz. He constructed and demon- 


in this case Rochelle salt—and has the primary crystal oscillator patent- 
Nicolson’s circuit is shown in Fig. 1.3. The crystal is effectively in a path 
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3 is in the direction of the smallest piezoelectric effect in Rochelle salt 
and contributes little to the action. If terminal one to the tapped coil is 
at the top of the coil, the circuit although employing a three electrode 
crystal connection, effectively reduces to B in which the crystal is in the feed- 
back path from plate to grid. On the other hand, if the tap is effectively 
at the bottom of the coil, the crystal is between grid and ground and feedback 
occurs through a distributed capacity from plate to grid. Both of these 
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Fig. 1.2—Ocean contour curves 


circuits B and C are widely used in oscillators of Pierce. Prof. G. W. 
Pierce published a circuit similar to circuit B, having a two electrode 
quartz crystal connected between grid and plate.’ 
In 1921, Professor Cady at Wesleyan University first showed* that quartz 
3 “Piezoelectric Crystal Resonators and Crystal Oscillators Applied to the Precision 
Calibration of Wave Meters,’’ G. W. Pierce, Amer. Acad. of Arts and Sciences, Oct 
1923, 81-106. 
‘“The Piezoelectric Resonator” W. G. Cady, Proc. ]. R. E. 1922, 10 83. 
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crystals could be used to control oscillators and that much more - 

oscillators could be obtained in this fashion. These were later applied to JB 282 
controlling the frequency of broadcasting stations and radio transmitters jy 9e adva 
general and about 1925 Mr. W. A. Marrison applied them to obtain a ver giver 
constant frequency and time standard, which is now used considerably by J elem 


the Bell System, by radio broadcasting systems, and by power companies not, 
The oscillators were subsequently improved by using crystals with sma cond 
temperature coefficients as described in Section V. At the present tim that 
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Fig. 1.3—Nicolson’s oscillator circuit of th 

crystal controlled oscillators are used very widely in radio military and alwa 

commercial applications. crys 

Another large use for quartz crystals is their use in providing very selectiv: as a 

filters. Probably the first use of a crystal to select a narrow frequenc\ dom 

range was made by Cady.‘ Using the very sharp maximum in curren! disc 

through a crystal at its resonant frequency, Cady proposed the use of such a chat 
crystal as a wave standard. This is equivalent to the use of a crysta 
as a tuned circuit. By incorporating a crystal in a three-winding trans- 
former and balancing out the static capacity of the crystal by an auxiliary 
condenser, W. A. Marrison® improved the selecting ability of a crystal used 


> Patent 1,994,658, filed June 7, 1927, issued March 19, 1935. 
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as a narrow band filter. At about the same time, L. Espenschied,® taking 
advantage of the knowledge of the equivalent electrical circuit of a crystal 
given previously by Van Dyke,’ showed how to combine other electrical! 
elements with crystals in ladder form to obtain band-pass filters. It was 
not, however, until the crystals were combined with auxiliary coils and 
condensers into the form of resistance compensated lattice type networks* 
that much progress was made in achieving the wide pass-band characteris- 
tics necessary for telephone and radio communication. Such filters have 


provided very selective devices which are able to separate one band of 


speech frequencies from another band different by only a small frequency 
percentage from the desired band. This property makes it possible to space 
hannels close together with only a small frequency separation up to a high 
frequency, and such filters have had a wide use in the high-frequency carrier 
systems, and in the coaxial system which transmits more than 480 conver- 
sions over one pair of conductors. In radio systems such filters have been 
used extensively in separating one side band from the other in single side- 
hand systems. 

In conclusion we can say, that the science of piezoelectricity was born 
about 62 years ago, lay dormant for nearly 40 years, but during the last 25 
years has advanced at such a rate that it can be regarded as one of the foun- 
dation stones of the whole communication art. 


1.3. THEORY OF PIEZOELECTRIC MATERIALS 


Piezoelectric crystals are of interest in communication circuits because 
they possess three properties. These properties are: (1) the piezoelectric 
effect provides a coupling between the electrical circuit and the mechanical 
properties of the crystal; (2) the internal dissipation of most crystals and 
particularly quartz crystals is very low, and the density and elastic constants 
of the crystals are very uniform, so that a crystal cut at a given orientation 
always has the same frequency constant; and (3), at specified orientations 
crystals can be cut which have advantageous mechanical properties such 
as a small change in frequency with a change in temperature, or a free- 
dom from secondary modes of motion. It is the purpose of this section to 
discuss the first property, the coupling between the electrical and me- 
chanical properties of the crystal. 

The piezoelectricity of quartz and other materials is due to the fact that 

* Patent 1,795,204, filed Jan. 3, 1927, issued August 8, 1933. 

*K.S. Van Dyke; Abstract 52, Phys. Rev. June 1925; Proc. I. R. E. June 1928. 

* See “Electrical Wave Filters Empioying Quartz Crystals as Elements,” W. P. Mason, 
B.S. T, J., Vol. XIII, p. 405, July 1934; “Resistance Compensated Band Pass Crystal! 
Filters for Unbalanced Circuits,” B. S. T. J., Vol. XVI, p. 423, Oct. 1937; “The Evolution 


of the Crystal Wave Filter,” O. E. Buckley, Jour. App. Phys., Oct. 1936; and Patents 
1,921,035; 1,967,249; 1,967,250; 1,969,571; 1,974,081; 2,045,991; 2,094,044. 
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a pressure which deforms the crystal lattice causes a separation of the cen- 


ters of gravity of the positive and negative charges thus generating a dipole 
moment (product of the value of the charges by their separation) in eac! 
molecule. How this separation can cause a coupling to an electrical circuit 
is illustrated by Fig. 1.4 which shows a crystal with metal electrodes norma 
to the direction of charge separation. If we short-circuit these electrodes 
and apply a stress which causes the centers of gravity of the charges t 
separate, free negative charges in the wire will be drawn toward the electroc 
in the direction of positive charge separation, and free positive charges i: 
the wire will be drawn to the electrode in the direction of negative charg: 
displacement until the crystal appears to be electrically neutral by any test 
conducted outside the crystal. When the stress is released the charges in 
the wire will flow back to their normal position. If, during the process 
we connect an oscillograph in the short-circuited wire, there will be a pulse 
of current in one direction when the stress is applied and a pulse in the oppo- 
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Fig. 1.4—Method for transforming mechanical energy into electrical energy in a crystal 


site direction when the stress is released. By putting a resistance in the 
connecting wire and applying a sinusoidal stress to the crystal, an alternating 
current will flow through the load and consequently mechanical power wi 
be changed into electrical power. Using the converse effect, a source of 
alternating voltage in the electrical circuit will produce an alternating stress 
in the crystal, and if this is working against a mechanical load, the electrica 
energy will be changed into mechanical energy. 

To apply this concept to quartz let us consider Fig. 1.5, which represents 
the approximate arrangement of molecules in a quartz molecule. Lord 
Kelvin’s explanation of the piezoelectricity of quartz is the following: 


“The diagram (Fig. 1.5A) shows a crystalline molecule surrounded by six 
nearest neighbors in a plane perpendicular to the optic axis of a quartz crystal 
Each silicon atom is represented by + (plus) and each oxygen double atom 
(minus). The constituents of each cluster must be supposed to be held together 
in stable equilibrium in viture of their chemical affinities. The different clusters 
or crystalline molecules, must be supposed to be relatively mobile before taking 














a 
tne 


infl 
cau 
the 


m¢ 
are 
at 
th 
in 
co 
eli 
me 
th 
to 


lit 
px 
in 
ro 






























Lord 


ving: 








QUARTZ CRYSTAL APPLICATIONS 185 


positions in the formation of a crystal. But we must suppose, or we may suppose, 
the mutual forces of attraction (or chemical affinity), between the silicon of one 
crystalline molecule and the oxygen of a neighboring crystalline molecule, to be 
influential in determining the orientation of each crystalline molecule, and in 
causing disturbance in the relative positions of the atoms of each molecule, when 
the crystal is strained by force applied from without 

“Imagine now each double atom of oxygen to be a small negatively electrified 
particle, and each atom of silicon to be a particle electrified with an equal quantity 
of positive electricity. Suppose now such pressures, positive and negative, to 
be applied to the surface of a portion of crystal as shall produce a simple elongation 
in the direction perpendicular to one of the three sets of rows. This strain is 
indicated by the arrow heads in Fig. 1.5A and is realized to an exaggerated extent 
in Fig. 1.5B. 





¥ 
Y A / Bee 
~ > 7s. DC 
+) t / 
¥ — ,)— 
+ — / 
© / / 
+)-C ' rw +)(— 
fa x ae f t yee ae (-F-) DHA 
+ % a + — +> 
/ 
> } 4 
/ j 
¥ + ¥ + 
+ + 
y fy 
(A) (B) 


Fig. 1.5—Kelvin’s model of quartz molecules 


“This second diagram shows all the atoms and the centers of all the crystalline 
molecules in the positions to which they are brought by the strain. Both diagrams 
are drawn on the supposition that the stiffness of the relative configuration of 
atoms of each molecule is slight enough to allow the mutual attractions between 
the positive atoms and the negative atoms of neighboring molecules to keep them 
in line through the centers of the molecules, as Fig. 1.5A shows for the undisturbed 
condition of the systems, and Fig. 1.5B for the system subjected to the supposed 
elongation. Hence two of the three diameters through atoms of each crystalline 
molecule are altered in direction, by the elongation, while the diameter through the 
third pair of atoms remains unchanged, as is clearly shown by Fig. 1.5B compared 
to Fig. 1.5A. 

“Remark, first that the rows of atoms, in lines through the centers of the crystal 
line molecules, perpendicular to the direction of the strain, are shifted to parallel 
positions with distances between the atoms in them unchanged. Hence the atoms 
in these rows contributed nothing to the electrical effect. But in parallel to these 
rows, on each side of the center of each molecule, we find two pairs of atoms whose 
distances are diminished. 
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“This produces an electrical effect which, for great distances from the mol: 
is calculated by the same formula as the magnetic effect of an infinitesima 
magnet whose magnetic moment is numerically equal to the product of the qua 
of electricity of a single atom into the sum of the diminutions of the two dista 
between the atoms of the two pairs under consideration. Hence, denoting | 
the number of crystalline molecules per unit bulk of the crystal; by 6 the radi 
the circle of each crystalline molecule; by g the quantity of electricity of each o 
six atoms or double atoms, whether positive or negative; by @ the change of dir 
tion of each of the two diameters through atoms which experience change of dire: 
tion; and by yu the electric moment developed per unit volume of the crystal, }y 
the strain which we have been considering and which is shown in Fig. 1.5B; we hay 


np = Nq4b0 cos 30° = 21/3 Nbq0 (1.1)” 


Kelvin’s model shows some of the symmetry properties of quartz. The 


axis marked X is the X or electrical axis of the crystal. The Z or optic axis 
is normal to the plane of the paper. The Y or mechanical axis is the axis 








(A) 


Fig. 1.6—Longitudinal and shear strains applied to a quartz molecule 


along which the stress is applied. It is obvious that if we rotate the dlirec- 
tion of the applied stress by 120°, a similar separation of charges at right 
angles to the stress will occur. There are then three electrical axes and three 
mechanical axes so that the optic axis can be regarded as an axis of threefold 
symmetry for the crystal. 

As can be shown from an extension of Kelvin’s model there are two other 
types of stresses that will produce a charge separation normal to the axis 
Suppose that we stress the crystal along the X or electrical axis as shown by 
Fig. 1.6A. Applying the same reasoning as before, we see that the apex 
molecules are separated farther apart without changing the separatio: 
between the other molecules. This results in a separation of the centers 
of gravity of the positive and negative charges, with the negative charges 
moving toward the left and the positive charges moving toward the right 
The separation is still along the electric axis, but is in the opposite direction 
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to that caused by a stress along the VY axes. A detailed analysis shows 
that the value of the electrical separation moment (dipole moment) for a 
stress along either axis is the same value but the sign is reversed. <A longi- 
tudinal stress then can only produce a charge moment along the Y or 
electrical axis which is the origin of the name electrical axis. 

If, however, we introduce a different kind of stress known asa shearing 
stress, a Separation of centers of charges can occur along the mechanical or 
Y axis of the crystal. A simple shear stress is one in which forces act normal 
to the direction of space separation rather than along it as shown, for exam 


ple, by the two opposed arrows normal to the mechanical axis in Fig. 1.6B. 


Such a shear does not occur in nature, but rather a pure shear which consists 
of two simple shears which are directed in such a way as to produce no 
rotation of the molecule as a whole about its axis. If we resolve these 
force components along directions 45° from the crystal axes, a pure shear is 
equivalent to an extensional stress along one 45° axis and a compressional! 
stress along the other 45° axis. Such a stress would cause the charges to be 
displaced from their normal position, as shown in the figure. This causes 
the center of positive charge to be displaced downward along the mechanical 
or Y axis of the crystal while the center of negative charge is displaced up- 
ward along the mechanical axis. 
These three relations can be written in the form 


P, = —dyX.+ dn, ; Py = 2dnX, (1.2) 


where P, is the polarization or charge per unit area developed on an electrode 
surface normal to the Y axis due to the applied longitudinal stresses Y, 
and Y,, while P, is the polarization normal to the Y axis caused by the shear- 
ing stress X,. dy is the piezoelectric constant and equations (1.2) show 
that the magnitudes of all these effects are closely related. In addition 
to these three major piezoelectric effects, quartz has two smaller effects 
which, since they are connected with the distribution of molecules in the VZ 
and XZ planes, cannot be demonstrated by the figures given previously. 
The complete piezoelectric relations are then 


i a = —dyX, + ayy, Dae duY, ’ P, _ dysZ, + 2d, X, tio3 


where Y, and Z, are respectively similar shearing stresses exerted in the YZ 
and ZX planes respectively. The best values for the di: and dy, constants 
are respectively 

g €.S.U.. 3 €.S.U 


; dy = 2.56 X 10 © (1.4) 
dyne dyne 


dy = —6.76 X 10 


as obtained by recent measurements for a number of XY cut and rotated X-cut 
crystals discussed in appendix A. 
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(Juartz is not the only type of crystal which is piezoelectric. Ip {a;; 
there are hundreds of crystals that exhibit this property. Whether 
crystal is piezoelectric or not and the relation between the stresses a; 
charge displacements depend on the symmetry of the crystal. Whenever 
there is a center of symmetry; that is, when the properties of the cryst, 
are the same in both directions along any line, no piezoelectric effect ¢, 
occur. This is illustrated by the simple arrangement of atoms shown } 
Fig. 1.7. It is obvious that no symmetrical application of forces can separat 
the center of gravity of the charges and hence such a crystal will not }, 
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Fig. 1.7—Crystal with a center of symmetry 





piezoelectric. Crystals can be classified into 32 possible classes on the basis 
of the symmetry they exhibit; and of these 32 classes, 20 are piezoelectri 
and 12 are not. As illustrated by the model for quartz, the response 1 
different types of force depends solely on the type of symmetry existing i! 
the crystal. 
1.4. ELEcTRICAL IMPEDANCE AND Low DIssIPATION IN CRYSTALS 

The first crystal used by Cady in controlling oscillators, was a crysta! cut 
with its major faces perpendicular to the X or electrical axis and with its 
length along the Y or mechanical axis. Referring to Fig. 1.5B, we see that a 
stretch along the Y axis will produce a charge displacement along the / 
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or X axis. Conversely, a voltage applied along the X axis will produce a 


iT 
} 
I 


sarge displacement and consequently a mechanical stress along the V axis 


( 
which will set up a longitudinal wave along the mechanical axis. As shown 
by Fig. 1.8, the type of motion resulting when the crystal is free to move on 
the ends is one in which the center is stationary and the ends move in and 
out. The crystal can then be clamped at its center or mounted from leads 
soldered to electrodes deposited on the surface. 

In using a crystal in an electrical circuit it is desirable to have an electrical 
equivalent circuit which will represent the electrical impedance as measured 
from the terminals of the crystal. Such a circuit’ is shown in Fig. 1.8. In 
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Fig. 1.8—Longitudinally vibrating crystal and electrical equivalent circuit 


Co 


this representation Co is the static capacity of the crystal which would be 
measured if the crystal were held from moving. C; is the stiffness of the 
crystal transformed into electrical terms through the piezoelectric effect of 
the crystal, while Z; is the effective mass of the crystal also transformed into 
electrical terms. The resonant frequency of the crystal is determined by 
the Young’s modulus and density of the bar according to the usual formula: 

fe = 5, Zs (1.5) 

2 p 


* Circuits of this type for representing the electrical impedance of a crystal were first 
derived by Van Dyke; see reference (7). The method of deriving them from Voigt’s 
equations is discussed in the appendix. 
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where Yo is the value of Young’s modulus along the bar, p the density, and 
f the length of the bar. 

A significant feature of the equivalent circuit is that there is always a 
definite ratio between Co and C; for a given crystal cut. This is really a 
measure of the ratio of electrical to mechanical energy stored in the crystal 
under an applied constant voltage. The reactance characteristic of the 
network is shown by Fig. 1.8 as a function of frequency. The reactance 
starts out as a negative reactance at low frequencies, becomes zero at the 
resonant frequency fx , becomes positive and very large at the anti-resonant 
frequency f, , then again becomes a negative reactance. Due to the high 
ratio of Co to C; existing in a crystal the separation between f, and fg be- 
comes very small. For example, for an AT crystal this ratio is around 200 
and the separation of f, from fz is only a quarter of a per cent in frequency. 
Since it can be shown that an oscillator will only oscillate on the positive 
reactance part of the crystal characteristic, the narrow separation between 
resonant and anti-resonant frequencies explains why a crystal can act as 
such a good stabilizer for an oscillator. As long as the crystal resonance 
itself does not change with temperature or other conditions, the very sharp 
reactance frequency characteristic will not allow the oscillator frequency 
to change much with a change in oscillator voltage, tube conditions, or any 
other changes which are likely to cause a change in frequency for a coil and 
condenser controlled oscillator. 

Strictly speaking, a resistance should be added in series with the jnduc- 
tance L, to represent the internal losses in the crystal, the loss of energy 


at the clamping points and the loss of energy due to setting up of air waves 


by the crystal motion. However, the value of this resistance and the 
amount of energy lost is very small in a crystal compared to what the losses 
are in purely electrical elements. A demonstration which shows this 
effect and shows that most of the losses of a well mounted longitudinally 
vibrating crystal are acoustic losses caused by setting up air waves in the 
vicinity of the crystal, can be made by using two oscillators, one a fixed 
oscillator and the other one controlled by a resonant circuit or a crystal. 
The fixed oscillator may be set at 99 kilocycles and the crystal oscillator 
controlled by a 100-ke crystal. The two will beat together giving the 1000- 
cycle note. When the battery is taken off the crystal] oscillator, it continues 
to oscillate till the energy built up in the crystal is dissipated in the interna! 
dissipation of the crystal. A good electrical circuit which has a ratio of 
reactance to resistance, or Q of the coil of 300 dies down almost instantane- 
ously. Fora crystal mounted in air it takes about half a second to become 
inaudible. This corresponds to a Q of 30,000 where Q is defined as the ratio 
of the reactance of the coil Z; of Fig. 1.8 to the resistance. For a crystal 
mounted in a vacuum a much higher (Q is obtained due to the elimination of 
the loss of energy by acoustic radiation. For such a crystal it takes eight 





secone 


S 1000 t 


As! 


tudin: 


alon 
othe 
mod 
the. 
Xo 
to e€ 
und 


QUARTZ CRYSTAL APPLICATIONS 191 


seconds to die down which corresponds to a Q of 330,000, which is about 


1000 times as great as that for a good electrical circuit. 
1.5. Mopes OF MoTION AND CRYSTAL ORIENTATION TO PRODUCE Low 
TEMPERATURE COEFFICIENT CRYSTALS 


As mentioned previously the first crystal cut used in oscillators was a longi- 
tudinal vibration along the Y or mechanical axis excited by a field applied 
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Fig. 1.9—Oriented quartz crystal cuts in relation to the natural crystal 


along the electrical or XY axis. This mode gives a good resonance free from 
other modes and a modification of it is now used in crystal filters. This 


modification, as shown by Fig. 1.9, (—18.5° filter crystal) consists in rotating 
the direction of the length by 18.5° from the Y or mechanical axis, about the 
X or electrical axis. As described previously”, the effect of this rotation is 
to eliminate the coupling between the desired longitudinal mode and the 
undesired face shear mode, thus simplifying the motion and eliminating an 


” “Electrical Wave Filters Employing Quartz Crystals as Elements,” W. P. Mason, 
B.S. T. J., Vol. XIII, p. 405 July 1934 or patent 2,173,589. 
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undesired resonance. However, to get a reasonably high frequency 
of a length vibrating type of cut requires too small a length to be practic: 

It was not long before crystal oscillators were controlled by thicknes 
vibrating crystals whose frequencies were determined by the thickness vf 
the crystals or by their smallest dimension. Referring to Fig. 1.6A, we «ff 


] 


that the same X cut type of crystal will generate a vibration along the ele 
trical or X axis when a field is applied along this axis. Since the thicknes 
dimension can be made very small, a high frequency is obtainable. Ho, 
ever, when the smallest dimension is used to control the frequency, a difi 
culty arises not present when the largest dimension is used to control t! 
frequency, namely, that harmonics and overtone modes of all the lower 
frequency types of motion produce frequencies near the frequency of th: 
thickness mode and it is difficult to pick out the desired mode. This wa; 
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Fig. 1.10—High frequency shear mode of motion 


especially true for the thickness vibrating X cut crystal and led to its 
abandonment in favor of Y cut crystals vibrating in shear. 

As seen from Fig. 1.6B, when a voltage is applied along the Y or mechan- 
ical axis, a shear vibration is produced which tends to change a square int 
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a rhombus. For a large plate in which the edge dimensions are large com- 
pared to the thickness, the motion occurs as shown by Fig. 1.10. For suc! 
a plate the motion is perpendicular to the thickness, which is the direction 
of transmission of the wave, and hence a shear wave is sometimes called a 
transverse wave. The frequency of such a wave can be shown to be 


_ 1. fem 
s=3,/% (16 


where ¢ is the thickness of the plate, ce. is the shear stiffness constant and p 
the density. The use of Y cut plates considerably improved the per- 
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formance of oscillators since the plates do not have as many secondary modes 
of motion as do the X plates. They have, however, one drawback. The 
frequency increases about 86 parts in a million for every degree Centigrade 


; increase in temperature. This requires regulating the temperature quite 


closely. 

In order to improve on the performance of the Y cut crystal, investiga- 
tions were made by Lack, Willard and Fair, Koga, Bechmann, Straube! 
and others" on how the properties of such crystals varied as the orientation 
angle of cutting blanks from the natural crystal was varied. As shown by 
Fig. 1.9, the crystals investigated all had one edge along the X or electrical 
axis with the normals making positive and negative angles with the Y axis. 
All of these crystals will have a component of field along the Y axis, which 
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Fig. 1.11—Frequency constant of oriented ¥ cut crystals 


will produce a shearing motion until the angles of cut approach 90 degrees 
from the-¥ axis. The smaller the angle A; the more strongly will the shear 
mode be driven. However, advantageous elastic relationships can be ob- 
tained by using oriented cuts. As shown by Fig. 1.11, Lack, Willard and 
Fair found that the frequency constant of a rotated crystal expressed in kilo- 
cycle millimeters varied with angle of cut and that there was a minimum fre- 
quency at +31 degrees and a maximum at —59 degrees. It was subse- 
quently pointed” out that these minimum and maximum points were signifi- 
cant angles in the elastic behavior of the crystal for they were the angles for 

1 “Some Improvements in Quartz Crystal Circuit Elements,” F. R. Lack, G. W. Wil 
lard, I. E. Fair—B. S. T. J., Vol. 13, pp. 453-463, July 1934; R. Bechmann—HF Techn. u 
El. Ak. 44, 145 (1934); I. Koga—Rep. of Rad. Res. i. Jap. 6, 1 (1934); J. Straubel, Z. tech 


Phystk., 35, 179, 1934. 
2 See patent 2,173,589. 
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which the high-frequency shear mode had zero coupling with the trou)!e. 
some low-frequency shear mode system of vibrations. Crystals cut at thes 
angles have a much cleaner frequency spectrum than Y cut crystals. Lack 
Willard, and Fair also found that the temperature coefficient of frequ: 


varied with angle as shown by Fig. 1.12. Starting from a high positive 
value for the Y¥ cut, the coefficient becomes zero at an angle of +35° — 15 
and —49°. The first angle cut is known as the AT cut and the second as 


the BT cut. Since the AT angle is nearer the Y cut, the piezo-electric con 
stant is larger and it is more strongly driven than the BT. On the other 
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Fig. 1.12—Temperature coefficients of oriented ¥ cut crystals 


hand, the BT has a higher frequency for the same thickness. Both crystals 
are near enough to the AC and BC cuts so that the systems of low-frequency 
shear modes are rather weakly driven. On the other hand, the shear mode 
of both crystals is rather strongly coupled to flexure modes of motion, as 
will be discussed by Mr. Sykes in a later chapter, and the crystal has to be 
exactly dimensioned in order that the flexure frequencies and other disturb- 
ing frequencies will not coincide with the desired shear mode. 

Other oriented shear crystals for lower frequency work are the CT and 


DT crystals investigated by Willard and Hight. They are related to the 47 


and BT crystals as shown by Fig. 1.13. The plate on the right shows the 
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motion of an AT plate. If we were to increase the thickness dimension until 
the plate was nearly square, the AT motion would correspond to a face shear 
mode which should be controlled by the same elastic constants as the AT 
motion. At the same time in order to drive the crystal efficiently we 
could decrease the width until it became the thickness. This procedure 
would be the same as cutting a crystal at right angles to the 47 and would 
suggest that by so doing we should obtain a low-frequency shear crystal 
with a low coefficient. Actually, Willard and Hight found that a crystal 
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Fig. 1.13—Relation of AT and DT low temperature coefficient crystals 


at —52° or 87° from the AT would give a low coefficient. This crystal was 
called the DT. Similarly, a crystal cut at +38° or 87° from the BT would 
also give a low coefficient and this has been called the CT. It can be shown 
that a component of the voltage applied along the mechanical axis will 
drive the shearing type of motion. The C7 is larger for the same frequency 
and more strongly driven than the DT. It is extensively used in controlling 
oscillators in the frequency range from 200 to 500 kilocycles. 

(Juite a few other crystal cuts have been discovered as shown by Fig. 1.9. 
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Another important cut is the GT,'"* which has a very constant frequency 


over a wide temperature range. As shown by Fig. 1.14, all zero temperatur, 
coefficient crystals are zero coefficient at one temperature only and usual| 
vary in a square law curve about this temperature. The GT crystal r pre. 
sented an attempt to get a crystal in which the frequency remained constant 
over a wide temperature range. As can be seen from the figure, when prop 
erly adjusted this aim is attained, for the frequency does not vary more tha; 
one part in a million over a 100-degree Centigrade range of temperatur 

This crystal makes use of the fact that a face shear vibration can ly 


resolved into two longitudinal vibrations coupled together. As shown by 














30 ) - 
z | | | | 
= 20 - | | | | } | |_| i 

A 

ae Cian 
fod 10 + aT + + + + . ay + ee = + 
a ie Fog | | 
a) hss ae) | | 
= 0 =r Se — + 5 — 
< eer \ ANa NN a ON i 
= ob rd et 
¥ 7 P X X | TY * | 











3 / [OOUGHNUT f on 7 OI 
50 at ‘ - oS a a ¥ . te-+ shames \- ——} 
ST Ti A STAY ye TN 
" / / Beer. \ \ | 
£ 5 | : — 1 ~ be \ if 4 + | a 
» s) / {CONG BAR, LENGTH i : 
Zz IXCALONG X AXIS. OS \ oN | \ 
S : te_T ST HARMONIC ’ } * T |* 
g sol |_/*~2Rb HARMONIC Ay Bi 8 RE lacs See. \ 
w / | | | | Er} 
i | ie Ea, 
ad: 20 30 40 $0 6 70 #60 #90 100 
TEMPERATURE IN DEGREES CENTIGRADE 
Fig. 1.14—Temperature frequency characteristics of a number of low temperat 


coefficient crystals 


Fig. 1.15, if we cut a crystal at an angle of 45 degrees from that of a shear 
vibrating crystal, an expansion occurs along one axis and a contraction 
along the other indicating that a face shear can be resolved into two longi 
tudinal modes that are coupled together. Now since it can be shown that 
all pure longitudinal modes for blanks cut in alli possible directions in a 
quartz crystal will have zero or negative temperature coefficients," it follows 
that if we have a shear vibrating crystal with a positive coefficient, that 

18 “A New Quartz Crystal Plate, Designated the G7, Which Produces a Very Constant 
Frequency Over A Wide Temperature Range,” W. P. Mason, Proc. 1. R. E., Vol., 28 pr 
220-223, May 1940 

‘4 This can be proved as discussed in the appendix by combining the Voigt expressions 
for the elastic relations in a crystal with the measured temperature coefficients of the six 
elastic constants. 
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coefficient must have been caused by the coupling between the two modes. 
As a result of this observation it follows that if we have a shear vibrating 
crystal with a positive temperature coefficient and cut another crystal at 
45 degrees from this crystal, the strong coupled mode which corresponds to 
the shear vibration will also have a positive temperature coefficient \s 
we grind down on the side, the two modes become farther apart in frequency 
and less closely coupled. Then, since they both will have a negative coeffi 
cient if separated far enough, it follows that for some ratio of axes, one of the 
modes will have a zero coefficient. This was tested out for a series of orien- 
tations near the CT and DT with the results shown in Fig. 1.16. Positive 
angle crystals had zero coefficients at ratios of axes varying from 1 to .855 
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Fig. 1.15—Relation between a face shear mode and two coupled longitudinal modes 
I 


depending on the angle while negative angle crystals had zero coefficients 
at ratios from .64 to 1.0. For positive angle crystals it was the higher fre- 
quency mode that was the stronger and could be given the zero coefficient, 
while for the negative angle crystals it was the lower frequency mode that 
was the stronger and corresponded to the face shear mode. 

Several of the positive angle crystals were measured over a temperature 
range with the results shown by Fig. 1.17. For angles above 51°-30’ the 
curvature was positive, while for angles below 51°-30’ the curvature was 
negative. Right at 51°-30’ the large square law curvature term disappeared 
and the frequency was constant to one part in a million over a 100-degree 
Centigrade range centered at 50°C. as shown by Fig. 1.18. Some further 
experiments showed that this flat range could be moved around a bit by 
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changing the angle of cut and the ratio of axes simultaneously. To go from 
—25°C. to +75°C. with a mid-range at 25°C., a crystal cut at 51°-7.5’ 
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Fig. 1.17—Temperature frequency relations for various angles for G type crystals 


with a ratio of axes of 0.859 is required. The GT crystal has been used 
quite extensively in frequency and time standards and in filters meeting 
rigid phase requirements. 
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Fig. 1.18—Temperature frequency characteristic for GT crystal 
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Fig. 1.19—Temperature coefficients of long thin rotated X cut crystals. Angle of 
rotation measured between length and Y axis. Dots are measured points. Solid line 
calculated from temperature coefficients evaluated in the appendix. 
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Two other cuts not previously described are shown also by Fig. 1.9. Thy 
are the MT low coefficient longitudinally vibrating crystal and the VP |o, 
coefficient flexurally vibrating crystal. Both of these are related to th, 
+5° X cut crystal of Fig. 1.9. As shown by Fig. 1.19 a long thin 5° ¥ cy: 
crystal is the best length direction for an X cut crystal to obtain a loy 
temperature coefficient. Figure 1.19 plots the temperature coefficients for 
long thin oriented X cut crystals, and this data is used in the appendix 
derive the temperature coefficients of the six elastic constants. However 
as the width of the crystal is increased the temperature coefficient becomes 
highly negative as shown by Fig. 1.20. 
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Fig. 1.20—Temperature coefficient of a +5° X cut crystal (g = 0°; 8 = 90°; y = 85 
as a function of the ratio of width to length. Ratio of thickness to length = 0.05. 


This change of coefficient occurs due to the fact that as the crystal width 
is increased, the face shear mode of motion becomes more strongly excited 
and contributes to the elastic constant. Then since the temperature coefii- 
cient of the shear elastic constant is highly negative for this orientation 
the temperature coefficient of the +5° X cut crystal becomes more highly 
negative as the width is increased. 

The MT longitudinally vibrating crystal employs a rotation of the plane 
of the crystal cut about the Y’ or length axis. The effect of this rotation is 
to change the temperature coefficient of the shear mode from highly nega- 
tive to nearly zero. The result is that the temperature coefficient becomes 
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very low and nearly independent of the width length ratio. The VT low 
coefficient flexurally vibrating crystal is similar to the M7 but requires a 
somewhat higher rotation about the Y’ axis to produce a low coefficient. 

The MT crystal has been used in narrow band filters such as pilot channel 
filters of the cable carrier system and in oscillators having frequencies be- 
tween 50 kilocycles and 100 kilocycles. The V7 flexure crystal is capable 
of producing frequencies as low as 4 kilocycles, and can be used to produce 
filters and control oscillators in the frequency range from 4 kilocycles to 50 
kilocyeles. Crystals of this type have been used with the Western Electric 
frequency modulation broadcast transmitter. Operating in the region 
of 5 kilocycles, they maintain the frequency of the transmitter to +.0025 
per cent without temperature regulation. These two crystals will be 
described in more detail in a subsequent paper. 


APPENDIX A 


VOIGT’S ELASTIC AND PIEZOELECTRIC RELATIONS AND 
THEIR APPLICATION TO THE DETERMINATION OF LOW 
TEMPERATURE COEFFICIENT CRYSTALS 


A.1 MATHEMATICAL EXPRESSIONS FOR PIEZOELECTRIC RELATIONS 


As mentioned in the historical introduction, Voigt formulated a mathe 
matical relation between the stresses, strains, polarizations, and electri 
fields existing in a crystal. For a general crystal devoid of symmetry these 
relations take the form 


sivX2 + ste Vy + sisZ. + sis V2 + sisZz 


—X, 
+ sisX, — du E, — du Fy — dn EF: 


Ey > -» Bee E = Bs 
—Vy = Sai Az + Soo } y Tt $23 Z; + $94 } z TT S25 bs 


tee 


Sze iy — daE, — dx ky, — dF: 


EK 


-} 
—i, = S31 », + S32 Ke -+ Sis Le + S34 1 535 y A 
+ 


s3eXy — dE, — dx Fk, — dE 
oat _ St Ne + S42 ry + S43 he + Si Fs + S45 vA 
+ sigX, — duEz — duF, — duk; 


“A New Broadcast Transmitter Circuit Design for Frequency Modulation,” J. | 
Morrison, Proc. J. R. E., Vol. 28, No. 10, Oct. 1940, pp. 444-449 
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soiXz + 552 Vy + 5532, + 554 Vs + 555 Zz 
+ sisX, — dE: — dx Ey — dyk 

—a, = 5X2 + seo Vy + sesZ. + ser Vs + se5Zz 
+ sésXy — dw Es — dx Ey — dy E 
P, = —dy Xz — du Vy — diyZ, — du Y; — dy Zs — dX, + «1 E, 
Py = —du Xz — de Vy — dyZ, — du VY, — duyZz — du X, +42 E, 


P, = — ds; aw —_ dse he — dz3 Z; — ds ie = ds, Zz —_ dy Xy ft K3 : 


where xz , yy , 2: are the three longitudinal strains, y, , 2, , x, the three shear. 
ing strains, X,, Y, , Z, the three longitudinal stresses Y, , Z, , X, the thre: 
shearing stresses; P, , Py , P, the x, y and s components of the polarization 
and E, , E, , E, the x, yand z components of the electric field. sf, ,--- , si, 
are the 36 elastic compliances. The superscript E is added to show that 
they must be measured when the field E is zero or the crystal plated and short 
circuited. As shown from section C of this appendix they can be measure: 
from the resonances of completely plated crystals. From the principle of 
conservation of energy it can be shown that there is the general relatio: 
between the elastic compliances 


E E 
Sij = Sji (A.2 


so that the greatest number of compliance moduli is 21. In equation (A.1 
the d;; are the piezoelectric constants measured by observing the propor- 
tionality between the strains and the applied fields in the absence of externa 
stresses. «x; are the “free’’ susceptibilities of the crystals in the three space 
directions measured in the absence of stress. The susceptibilities are related 
to the “free” dielectric constants K? by the equation 


K? =1+ 4x? (A3 


In addition to these equations we have also that the charge per unit area () 
on the surface is related to the field and polarization by 
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These equations hold for the most general type of crystal. In addition 
Voigt showed that if there was any symmetry existing in the crystal, a num- 
her of the constants were zero and certain relations existed between other 
constants. For example quartz has a trigonal symmetry about the Z or 
opti axis, and three digonal axes of symmetry (the three XY or electrical 
axes) about which it is necessary to turn through an angle of 180° before 
the original pattern is restored. Voigt showed that by expressing the rela- 
tions (A.1) in terms of rotated axes and imposing the symmetry condi- 
tions, the following relations existed between the elastic and _piezoelectri 
coefficients 
5 = S16 = 
st} 53: S13; 534 i — sts; 
S56 - 2st4; See = 
dy = diy = dy = dx dog = dau = dx, = dy 
1 short 
asure = dy = dy = d35 = dy = 0 
ple of dy —dy; dog = —dia; do = —2d 


lation 


Hence the relations between the stresses, strains, polarizations and fields for 


( ) . ° 
(A.2 quartz reduce to the simpler forms 


Ew |, E y- , Buy . 
(Al —iz = CD, oF a i Sjo)} y 51322 _ Sis) ae dyk, 
ropor- ] 


& r E , E , , . 
= SX, + Sit] y 7 $13Z SuVY, + dunk, 


terna 
“ie Ze = 53X22 + Siz3Vy + 5332, 
Plated 
= s4X,—- siuVy +- s0.¥, ~— dike 
= siuZe + steXy + duky 
2suZ-+ 2(sn — sty) X, + 2d),F, 
_ E,Ki 
4r 


. >F 
i “es ' + duZz + 2duX, 


_ EK; 
4r 


— dX, + du Vy — duY, 


The superscripts have been left off the constants s13 , 533 and K3 since it 
will be shown that their values are not affected by the way in which they 
are measured. 
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Equation (A.5) is not the only way of relating the elastic and ele 
tric quantities. For example if we substitute the values of the fields of thy, 
last three equations of (A.6) in the first six equations, we can write 


le == shi X, + sioVy + $132, + siuV, — guQz 
—yy = St2X2 + stiVy + 5132, — staVs + guQz 
—2, = 3X2 + Sis Vy + 5332, 

—ye = StuXz — StaVy + suVe — guQz 


Qz Qy 
—22 = SaaZz + 2sigXy + BuQy 


—x, = 2sfsZ, + 2(st) — sie)Xy + 2en0, \7 
6 4r : ‘ " 
E; = — 0; + guXz— guy + gu; 
Ky 
: 4x , 
Ey, 7 F 7 guZz = 2guXy 
Ky 
4r 
Ek. ==—Q: 
K, « 
where 
. And} , , 4ndi ,  4ndy d, 
. 11, Oo . i a 11 14 
a1 = oy K? } Siz = Siz + K? ; a4 @ ae kK! 
Q E dads 4rd, 4rd, 
6. * Su > 5 a= 5 i 
1 Ky Ky 


The superscript Q is added to show that these are the elastic compliances 
that will be measured when the free charge on the surface is zero. These 
elastic constants are the ones measured when an unplated crystal is put in 
an airgap holder with a large air-gap since then no charge can flow to the 
surface of the crystal. The difference between the zero field and zero charge 
elastic constants for quartz is less than 1 per cent. For rochelle salt, how- 
ever, they may differ by a factor of 4. For rochelle salt the principal piezo- 
electric constant dy and the “free” dielectric constant Ki vary widely i 

value and phase angle with variations in temperature and frequency, whereas 
the piezoelectric constant gi4 which is proportional to the ratio of these two 
is nearly a constant for all frequencies and temperatures, so that the formu- 
lation of equation (A.7) is more advantageous than that of equation (A.6 

For quartz, however, both forms are reasonably constant. Furthermore 
the elastic constants of equation (A.6) are those for a plated crystal which 
are usually the ones of interest for a crystal employed in an oscillator or filter 
Hence this formulation has been used in this appendix. 
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‘used later in this appendix, they are given below. 
Fobtained directly from equations (A.6) by solving them simultaneously to 
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Both the formulation of (A.6) and (A.7) can be expressed in terms of 
the strains rather than the stresses. Since these are useful forms and are 


Equations (A.8) are 


replace the strain by the stress, while equations (A.10) are obtained in the 


same way from 


equations (A.7). 


i ae oF — E 9p 
—X, = cn, + Ci2Vy TT C1322 TF C14Vz2 — enk, 


Y, = Ye a oer : ae E 
— Vy = (12Xz2 T CiiVy T C1332 — Craze T Crile 


—Z, = C13Xz + C13 Vy + C3322 





, E E E . 
—Y, = CuaxX: — CisVy T C44Vz — euk, 
r E ' E 7 
—Z, = €442r T Ci4Xy + euky 
cf ct 
‘3 E 11 — C12 . 
—Xy = 614% + ( ; )x + en Ey (A.8 
BE. FE, Ks 
QO: = in + P,= in + ¢1%3 — C11 Vy -- C14 Vz 
E E Ki 
WP=2—+P=— — 4%: — eux 
Cy dr 7 dr " 
0, = a E.Ks 
dr 4r 


where the relations for the elastic constants are 


E E E 
7B _ S33 S44. > -& _ 533 S44 ; = Sis 
a <9 “tig =: —_ } cE == ; 

a B a B a 

E E 
a. mn —S14 po Sir + Siz , a Sit — Si2 
ie 3 ’ “a > > 

B ’ a B 

E 

C6 ps — OR v Ses(St1 + Si2) si 
ua = ) = Sit TT Si2) — £513 
66 D 28 33(Si1 12) 13 
E E E E2 
B = Saa(Sit Si2) — 2sy4 
Conversely we can also write the useful relation 
7 E E E 
C33 Cas E C33 C44 
ds = — * a io. = = . 
#911 ? ,? #912 7 ? 
a B a B 
E E LE 
2. Sou i _ ou + C2 
$3 = — ys $14 = = $33 = r = 
B at 
Prd 
E Cur — C12. -. el BE, 444. 
S44 = 7 ; See = 2(S11 — Siz) = —; ; 
B 8 

E E 2 E k f R2 

a= €33(C11 -+ a3) = 2C13 ; B’ = CqyalCiu1 — C12) — 2C14 
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For the piezoelectric constants 
; E I a ‘ —_ _& - 
én = du(ein — ci) + ducis 5 14 = Qdirceig + dyacyy : 
and conversely 
E 
l 


E : ‘eo ae eae se E 
—dyy = (Sis — Siz) + eraSi4 5 —dyy = leuSi4 + er4544 


rr . . -C . . 

rhe dielectric constant A; denotes the clamped dielectric consta: 

the constant measured when the crystal is free from strain. This is relat, 
° . » -F ° 

to the free dielectric constant K, by the equation 


Ki = Ki - 4m[diseiusg + 2dureni]. 19 


In a similar way if we solve equations (A.7) simultaneously, for the stresses 
in terms of the strains, we have 


—X, = chxs + choyy + Crste + Clays — fz ; 
—Vy = chore + chivy + ct. — Clave + furOz ; 
—Zs = Ci3Xz + CisVy + C332: 5 

—V, = chute — chiyy + chiye — fuQz ; 


5 .Q Q P 
—Ze = Caz + Ciaty + fiQy ; 


ae. chi — Ch ’ = 
—Xy = cits + . xy + furQ,y ; (At 


4r O 


i 2 — Sute + futy — Sure; 


, 4 
E, = Pa + fist. + furry } 
1 


Ro 6. 
Rs 


where the c@ constants are related to the s@ constants as in equation (A.8 
The piezoelectric relations are 


2 .Q we < rs Q a 
fr = gu(cn — C12) + 21aCi4 5 fis = 2guer4 + 814044 5 


or conversely 


—" fu(sii — Sir) + eusis ; =f = 2fur sta + fis Si 5 
while the dielectric relation between the free and clamped crystal 
4 dr 
= — (gisfis + 2gu fu). (A.11 


n 
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Equations (A.10) might also have been obtained dire tly from equations 
(A.8) by substituting the charges from the last three equations in terms of 
the fields. This substitution yields the additional relations 


° 
4re}; 4rre}; 
e .. ll, 7 a wei) | Q E 
Cn = C1 + as os = os = 9 Cig = C13; 
Ky Ky 
iT 
: Q E drei C14 | Q E. 
relate Cis = Ca + 7 C33 = C33; 
i Ki 
k E 8ren? (A.12) 
2 .Q Q ws Or -C 
a. 4reiy @ C11 — Ci2 K, 
Cia = Coa + a Cea = = : 
Stresses RE; 2 ? 
4 dor 
fu Fa ass fia aeliwsr ge ie 
Ky Ky 


(A.2). VALUES OF THE ELASTIC AND PIEZOELECTRIC CONSTANTS 


The first and one of the best determinations of the elastic constants of 
quartz was made by Voigt. Using static deformations of unplated crystals 
he determined the elastic constants to be 


én = 85.1 X 10'° dynes/cm’; Ce = 6.95 XK 10"; 
(All (ig = 14.1 X 10"; C4 = 16.8 X 10!° 
C33 = 105.3 X 10"; ty = Sak 10°: (A.13) 


cos = (* > 2} = 39.1 x 10” 


From these the moduli of compliance can be calculated and are 


Su = 129.8 K 107" cm?/dyne; Sie = —16.6 X 10-4; 


53 = —15.2 x 10 sin Siu = — 43.1 x 10-14: 

aa (A.14) 
$33 = 99.0 X 10-**. Su = 200.5 X 10714; 

See = 2(Su — Sy) = 292.8 K 10-4. 


Whether these are zero field or zero charge constants is not known. If 
they were measured in a room with high humidity, the polarization produced 
; by strain would soon be annulled by a current flow through the leakage re- 
sistance of the adsorbed moisture, and the constants would be ct; or st; . 
On the other hand if the displacements were measured in a very dry room, 
(AM the leakage resistance is very small and it may take hours to annul the polari- 
zation through a leakage current flow. In that case the constants measured 
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would be c?; or s?; . In any case the difference was probably less thay thy 
accuracy of measurement. 

Later measurements by Perrier and Mandrot for two of the constants 
and 533 give the values 


$a = 1233 X 10™: Ss3s = 97 X 10-4 \.15 


\ 


By using the measured resonance frequencies of known modes of motio 
the uncertainty of the type of elastic constant can be removed, for th, 
alternations occur so fast that the leakage resistance has little effect. Ii 
crystal is lightly plated, it is shown in the next section that the resonant 
frequency of a length vibrating bar will be determined by the zero field 
elastic constants st; ._ On the other hand if an unplated crystal is measure: 
in an air gap holder with a large air gap it has been shown that’ the fr 
quency measured will be determined by the zero charge elastic constants 
si; or c¢;. A careful measurement of the elastic constants of quartz has 
recently been made by Atanasoff and Hart®. Using thickness modes for 


i 

1 The resonances of length vibrating crystals have been discussed by Cady, ‘‘The Piez 
electric Resonator and The Effect of Electrode Spacing on Frequency,” Physics, Vol. 7 
No. 7, July 1936, pages 237-259; and by the writer, “Dynamic Measurement of The ( 
stants of Rochelle Salt,” Phys. Rev., Vol. 55, pages 775-789, April 15, 1939; whilk 
resonances of thickness vibrating crystals have been discussed by Cady (above paper) ar 
Lawson “The Vibration of Piezoelectric Plates,’ Phys. Rev., Vol. 62, July 1, 1942 
71-76. Foralength vibrating crystal Cady shows that the resonant frequency for 
gap (plated crystal) is controlled by the constant 1/s?,. For a crystal with a larg: 
gap, the frequency is controlled by the constant. 


E 2/;KF Bb? _ 14/,@ 
1/sy, + 4nd,i/K, Sy, = I/sy,. 
Starting with equations of the form (A.10), the writer showed that the frequency of a bar 


an air gap holder would be controlled by the constant 1/s/,, while the frequency of a plat 
crystal is determined by 


For a thickness vibrating crystal for which the field is applied in the direction of wa 
propagation, Cady and Lawson find that the resonant frequency is controlled by the elast 
constant 


4re,;" 8 


~ io 
r rear (P a 


where D is the total separation between electrodes and ¢ the thickness of the crysta 
When the separation is infinite, the controlling elastic constant is cf, + 4me1,?/K{ whi 
from equation (A.12) is c’,. When the air gapis zero or D = 1, the controlling constant is 


> 4re,; 8 
an + Ke (: *) 
1 


which, for all practical purposes, can be taken as ¢;7 for quartz. 


?“Dynamical Determination of the Elastic Constants and their Temperature Coet! 
cients for Quartz,” Phys. Rev., Vol. 59, No. 1 (85-96), Jan. 1, 1941, 
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relatively thick pieces of quartz, and determining the asymptotic value for 
high order harmonics, they obtained the elastic constants 


_~ 


87.55 X 10'° dynes/cm*; Cie = 6.07 K 10! C3 = 13.3 & 10°: 
—¢4 = 17.25 X 10"; ‘93 = 106.8 XK 10!°; cu = 57.19 XK 10°. 
(A.16) 


In addition they came to the conclusion that css had a value of 18.4 & 10", 
which was different from the value of c4 as required by theory. Their 
measurements were made with high harmonics in air gap holders so that the 
values measured should determine the c%; constant. To explain the dis- 
crepancy found, Lawson® has suggested that the cy, constants 


Q KE 


Cijy = Cij TT 4rrey,e) K\ AAZ 


do not obey the same symmetry relations as the « * constants. This sugges- 
tion does not seem to be borne out by equations (A.10), from which the sym- 
metry relations of the c;, constants can be determined. If we start with a 
generalized form of these equations applicable to any crystal and apply the 
symmetry relations for quartz, we find that it is still necessary to satisfy the 
symmetry relations between the constants found previously and in particular 
Coe = C14 (A.18) 
In order to investigate this matter further, and to obtain more reliable 
values of the elastic constants, an analysis has been made of a number of 
measurements previously obtained for oriented crystals. In particular two 
families of oriented crystals were investigated. One family was a set of 
oriented X cut crystals vibrating longitudinally. They were cut with their 
major faces normal to the X axis and with their lengths at angles A» of from 
+43° to —79° with respect to the Y or mechanical axis. They were 
oriented similarly to the +5° and —18.5° filter crystals shown by Fig. 1.9. 
When these crystals are 7 to 10 times as long as they are wide or thick it 
has been shown previously* that their length resonances are determined 

very accurately by the equation 
he = 7 (A.19) 

eal 


gE 
pS22° 


where ¢, is the length of the crystal, p the density and sso the inverse of 


Young’s Modulus along the length for a plated crystal. This is related to 
the angle of cut Az by the equation 


E E 4 E - 4 Ez 3 . 
S92 = $i, COS Ao + S33 Sin Ae + 2514 cos As sin Ag 


(A.20) 
By <9 2 
a (2513 + s4,) sin’ Ae cos Ae 
\. W. Lawson, Phys. Rev., 59, 838 (1941). 
‘“Electrical Wave Filters Employing Quartz Crystals as Elements,’’ W. P. Mason, 
B.S.T.J., Vol. XIII, pp. 405-452, July 1934. See Figs. 25, 31 and 32. 
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Since the resonant frequency of the plated crystal was measured, it was the 
zero potential elastic constant that was determined. These crystals were 
lightly plated with aluminum and it had been previously shown that the 
added plating would affect the frequency by considerably less than 0.1 per 
cent. The crystal orientations, their dimensions, the frequency constants 
and the values of s%- are shown by Table I. 

These measured variations satisfy equation (A.20) for the variation of 


E . of 
22’ With angle very well if we take 


5 


127.9 x 10" cm’, dyne; 533 = 95.6 K 10 


’ 
—44.6 x 10°"; 
- 2513 = 175.8 X 10°". 
TABLE I 


Dimension, mm 
Resonant | Frequency E 
Frequency | Constant Value of s54, 
re yet Thick- 25°C KC cms 22 
Width 


ness 


Angle of 


Length 


24.03 
20.00 
19.99 
19.98 
20.02 
20.03 
19.97 | 
20.03 | 
19.96 | 
20.02 
19.97 

- 4.05 | 19.95 | 
+11.8 19.11 
+18. 20.02 | 
+42. 20.00 | 


.90 | .502 | 130,700 | 314.1 .6 X 107 cm?/dyne 
50 -902 | 127,710 | 255.4 5 
97 .505 128,390 256.8 A 
95 .500 | 128,590 257 .( 
92 500 | 132,130 264.: 
02 | .502 134,050 269. 

99 .502 | 135,240 | 270.: 
03 508 | 135,890 | 272.( 
98 506 | 135,920 272.1 
.98 505 136,890 274.0 
00 .505 | 138,400 276.5 
97 | 510 139,900 279.0 
O01 | .500 54, 295.4 
95 500 bP 311.1 
95 .500 ta 349.5 
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This gives three of the constants directly, and a relation between two more. 
To obtain the remaining constants and to test out the hypothesis that there 
are seven elastic constants rather than six, use has been made of measure- 
ments made for thickness vibrating shear crystals obtained by rotating one 
edge about the X axis. These are the AT and BT series shown by Fig. 1.9. 
As shown by a former paper®, the frequency of such crystals depends on 
the edge dimensions as well as the thickness dimensions. However, as the 
edge dimensions become large compared to the thickness dimension the 
principal frequency approaches an asymptotic value which is taken as that 
for the infinite plate. For the AT, BT and Y cut crystals these asymptoti: 
values have been determined to have the values shown by Table II. 


5 “Tow Temperature Coeflicient Quartz Crystals,” B.S. T. J., Vol. XTX, pp. 74-93, 
Jan. 1940. See Fig. 5. 





It we n 


liffers | 


The de 


The ot 


Tot 
can mé 


these € 


F Substi 


Comps 
they a 
not in 
six. | 
by Ta 

Mei 
electri 
stande 


 detern 


r of a | 
EB 
Soo ar 


QUARTZ CRYSTAL APPLICATIONS 211 


if we make the assumption that there are seven elastic constants and ce, 


Al 


\iffers from ¢i4 , the frequency of this series of crystals will be 
1 ay E FE 2 F E 
foe ~66" where Cee: = Ces COS Ao + C44 Sin Ao — Csgsin2A4. (A.22) 
)f 
- p 
lhe determination for the Y cut gives directly 
m] z 10 
C66 = 40.5 X 10° dynes per square cm. 
[he other two cuts give the values 


css = 18.2% 10°; ck, = 58.65 & 10” (A.24 


To test out the hypothesis that css differs from cts or sf, from 2s*, we 
can make use of equation (A.8) writing cg. in place of c{;. Then solving 
these equations simultaneously we find 

Cee —s — Cis 


> 356 = —> E2 


is BS E2 So ae ; 
(Cs4C66 — Coe ) (C4aCo6 — Coe ) 


TABLE II 


Asymy totic Fre 
Crystal Angle of Cut Az quency constant 
KC mms 


AT +35° 15’ 1661.5 29.39 & 10" dynes/cm? 
Y Cut 0 1954 40.50 
BT — 49° 2549 68 . 86 


Substituting in the values from (A.23) and (A.24) we find 
” 14 2; E 14 
197.8 X 10 ~ cm'/dyne; Sse = —89.0 X 10°: 
{ ) 
P E E > ind \.26 
= 2(Si1 —- $}2) = 286.5 XK 10 


Comparing the value of sj with 2s}, given in equation (A.21) we see that 
they are equal within the experimental error, so that these measurements do 
not indicate that there are seven elastic constants but only the customary 
six. Using these values all the elastic constants can be evaluated as shown 
by Table ITI. 

Measurements have also been made to determine accurately the piezo- 
electric constants. This was done by using the ratios of capacities of two 
standard rotated X cut crystals for which these ratios have been accurately 
determined. As shown by section C of this appendix, the ratio of capacities 
rof a crystal is related to the piezoelectric constant dx , the elastic constant 
sy and the free dielectric constant Kt by the equation 


, a r(i—k er 
r = ratio of capacities = — ( . ) (A.27) 
8 k? 
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where k the electromechanical coupling is given by 


, 4r 
acs Ty eee 


The two crystal cuts and their constants are given in Table IV, 
. . . , 
Only the numerical value and not the sign are determined for d> . 


TABLE III 
Elastic Compliance Moduli Elastic Stiffness Moduli 


= 127.9 X 10- cm?/dyne cf, = 86.05 X 10! dynes/cm 
sf, = —15.35 cf, = 5.05 

= —11.0 C13 = 10.45 

= —44.6 cf, = 18.25 

= 95.6 C33 = 107.1 

= 197.8 cf, = 58.65 
hh — 4: 0.5 


2(sf, — sf) = 286.5 ck, = 2 


TaBLe IV 
Angle of Cut, Az cunt. Value of a Value of xP Value of qd), 
—18.5° X cut 137 144.5 
0° X cut 125 127.9 
TABLE V 


Piezoelectric Value in cgs electrostatic Piezoelectric Value in cgs electrostat 
constant units constant 


dy .76 xX 10-8 Riu 
dis .56 x 10-8 Ria 
ei ‘ 01 x 104 fir 
C14 .97 x 104 tis 


The variation of di: as a position of angle has been shown to be® 
diy = —4[du(1 + cos 242) + dy sin 2A9] 
The two values of d,2 of table IV are satisfied by 


dy, = —6.76 X 10: dy = +2.56 X 10-8 (A.30 


® See “Electrical Wave Filters Employing Quartz Crystals as Elements,’ W. P. Mason, 


B.S. T. J., Vol. XIU, 405 (July 1934). 
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From these values and the elastic constants of Table III we can calculate 
all the different forms of the piezoelectric constants. These are given in 
Table V. 

(A.3). DERIVATION OF EQUIVALENT CrRCUIT OF CRYSTAI 


The electrical impedance and electrical equivalent circuit for a fully 


plated crystal can be derived from the piezoelectric relations of equation 


(4.6) taken together with Newton’s law of motion 

. dé 

F, = ma = (p dx dydz) — 
df? 
where m is mass of an elementary volume dx dy dz, a the acceleration, and 
¢ is the displacement of the element in the y direction. If we consider a 
long thin X cut crystal with its length in the y direction, the piezoelectric 
relations of interest are 


gE , gE r E , + 
—Yy = 82 Xz + Si Vy + 51322 — Sin V2 + du Ex; 
" >F ) 
E,K? : : : (A.32) 
—_ dX, +d,Y, meats dis Vz. 
4 
For a long thin crystal with its long dimension in the Y direction we can set 


X,=Z,=¥Y,=0 (A.33) 


(A.31) 


0. = 


This follows since the crystal is free from external forces, and hence these 
stresses on the edges of the crystal must be zero. On account of the small x 
and z dimensions, the rate of change of these ‘stresses with x or z will have 
to be high in order that the stresses shall differ appreciably from zero, and 
there are no mechanical strains causing a high stress gradient. Then for a 
long thin bar the piezoelectric equations can be written 
— yy, = sh Vy + dy, E,; 
E,K? , (A.34) 
de + dy, } ye 


Let us next consider a small cross section of the crystal with a dimension 


0, = 


dy along the crystal length. The total force on the section is a resultant of 
the difference in stresses on the two faces or equal to 


; ay 
fol ¥y, — Val = — bales dy = Fy (A.35) 


where Y, the stress is considered as a compressional force acting on the faces 
of the element. By Newton’s law of motion (A.31) we have 


oY, dt aY, dé 
— ly , = ly ; = — _ (A.3¢ ) 
ly t, dy ay lw f, dy p r or ay P op ) 
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For a completely plated crystal such as we are considering, the potent 
gradient E, will be independent of the y direction, since any charge distri} 
tion will be equalized with the speed of light which is much higher tha 
the speed of sound in the crystal. Then equation (A.34) when differentiat; 


by y becomes 


. 2 ” 
OY, “ og — OVy A 3" 
_— = - Sit ‘ (ALS 
oy oy" ov 
Introducing equation (A.36), the equation of motion for a plated crysta 
becomes 
2 2 
vg s @é ' 
ae SE SRD aes (A.38 
oy" dt? 


For simple harmonic motion the variation of — with time can be written i 
the usual form 


wt > 
g Ping ge” , (A.39 


so that for simple harmonic motion equation (A.38) becomes 


dé 2g d w 
> — w Sit pc = $ —--—t=0 (A40 
dy" dy’ v 
where v the velocity of sound in the plated crystal is given by the formula 
2 1 ‘ed 
v,=.—. (A.41 
pSi1 


A solution of equation (A.40) with two arbitrary boundary conditions is 


" w tt ; 
§ = Acos—- y+ Bsin Z (A.42 
e* v 
To determine the constants A and B, use is made of equation (A.34 
Differentiating (A.42) 


0é ) — | wy Pp . . 
—-—-=-y= ~| Asin— y — Bcos— | = Sh Y,+dnk, (A343 
oy ° v o v 


When y = Oand y = ¢ the bar length 


Ys = Fo and Y, = Y,, (A.44 


provided the crystal is driving a load. For most electrical cases the on)\ 
load driven is an air load and this is usually very small so that it is customary 
to set Y,, = Y,, = 0. Under these conditions 


—" B= dF, and | 4 sin wot — Boos =4 = dyE,. (A.45 
v v v v 
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Solving these equations for A and B and substituting in (A.43) we have 


wl . wy wy — : 
5, Sin + cos | = $1 7, + dy FE, 
£U v v 


— wly — €/2) (A 46) 
7 hel cos > | 


or Y, — E 1 = . 
Si | wl | 
cos 
2v 


The electrical impedance measured at the terminals of a plated crystal is 
then determined by substituting the value of Y, in the last of equations 
(A.34) and integrating the charge Y over the whole surface. The current 


| 
~ 
© 

| 


= di; Ez tan 


into the crystal is then 


-P 2 F : 
Ky diy cos w(y — 4 2) 
= -> 1— = dy 
4r Si1 wl 
cos 


{ 
i= joQ = jole | Es 


= juE, tot) — —-S11 -—— , 
jul z | (A.47) 


2v 


ft 
es } 9 tan Sf] 
| ts 


>LC >F 4dr diy P P . ° . 

where Ki = Ki — ~~ 3S called the longitudinally clamped dielectric 
Sil 

constant, i.e. the dielectric constant that would be measured if we suppress 

the longitudinal strain along the y axis but not the other strains. The 


admittance of the crystal then is 


| tan “| 
P ° ‘ >LC 9 la a 
Bap 2 ad jut, €| Ky 4 du 20 : (A.48) 
e” Ze Site" 2 a | 
2v 


This consists of two terms which represent parallel branches in the equiva- 
lent circuit. One of these is the capacitance 


_ ty (Ky° t, (Ky ° 


Co = txt, cgs units = tnt,9 X 108 farads (A.49) 





216 BELL SYSTEM TECHNICAL JOURNAL 


The other branch contains the impedance 


wl E wl 


° ll 
—jt,| sh 2s “ag XP KO 


‘e 20 
> cgs units = ohms 
wl, f dy wl 2 Ww 
tan - wl, f dy tan 


dy 


This branch will have a zero impedance or will resonate when the tangent js 
infinite or when 
2rfel T fi a | 1 
=-—-oOrf = —= (A.51 
2v 2 2 2 k ‘a 
UV pst 


Hence for a fully plated crystal it is the zero field elastic constant that 
determines the resonant frequency. 

Near this resonant frequency, the impedance of equation (A.50) can be 
represented by a series capacitance and inductance having the values 


ae fof 8 ‘.” I = ps*, tt, x 9 b 4 10" 


Cy > ; , <1 » 
“ws x9 x 10" et. du” 


(A.52 
Taking the ratio between Co and C; we have 


‘' 4r a)" 


Oe =( : " _® si Ky 
C; 8 dr 4." 8 4a d;? 


>F EB 
Ky Sit / 


where & the coefficient of electromechanical coupling is equal to 


4r P 
b= dng/ os Bs (A.54) 
| Ki si 


These values are used in equations (A.27) and (A.28) to evaluate the piezo- 
electric constants of quartz. 


A.4. Use oF Vorct’s RELATIONS IN LOCATING REGIONS OF LoW TEMPERA 
TURE COEFFICIENT CRYSTALS FOR SIMPLE MODES OF MOTION 


In Section 1.5 of the text, the statement is made that all longitudinal] 
vibrating crystals of quartz have a zero or negative temperature coefficients 
This can be proved from Voigt’s relations for quartz and a knowledge oi 
the temperature coefficients of the six elastic constants of quartz. Since 
the same method can be used to locate the regions of low temperature coefli- 
cient for other simple modes of motion a short discussion of the method is 
given here. 

The Voigt relations given in equation (A.6) give the values of the piezo- 
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electric and elastic constants for crystals with their three edge dimensions 
along the three crystallographic axes. Most low-coefficient crystals, how- 
ever, are oriented crystals with one or more of their edges lying along 
directions not parallel to the crystallographic axes. The theory of elasticity, 
however, provides methods for calculating the values of the constants for 
rotated axes. If the rotated axes X’, Y’, Z’ are related to the crystallo- 
graphic axes X, Y, and Z by the relation 


XV¥Z 


fm m (A.55) 
" fy M2 Nz 


| (3 ms Ms 


where ¢; , «++ , #3 are the direction cosines between the axes indicated, the 
theory of elasticity provided relations between the stresses of the rotated 
axes and the stresses of the crystallographic axes, between the strains of 
the rotated axes and the strains of the crystallographic axes, and between 
the field, polarizations, or charges of the rotated axes and the same quantities 
for the crystallographic axes. Then if we express’ the relation between the 
stress, strain and fields for the rotated axes, the elastic and piezoelectric 
constants are determined. 

Two shorthand methods are also available for calculating the constants 
of rotated crystals. One method” is the matrix method which is based 
upon the fact that relations in (A.6) can be expressed in a matrix equation 


—e=s*X +dE (A.56) 


where ¢ are the strain components, X the stress components, s” the elastic 
compliance matrix, d the piezoelectric matrix and E the field components. 
By applying the rules of matrix multiplication the s and d matrices can be 
transformed to rotated axes having the direction cosines of equation (A.57) 
with respect to the crystallographic axes. The other method is the method 
of tensor analysis. Equations (A.6) can be expressed in the form” 


E . ‘ ' _ 
—— SijapX ap a dijpEy (A.57) 


where ¢;; is the second rank strain tensor, Xag the second rank stress tensor, 
Sija3 the fourth rank compliance tensor, F, the field vector, and d;;, the third 
rank piezoelectric tensor. By employing the geometric rules for tensor 


7 This method of determining the constants for rotated axes is discussed in a former 
paper ““Dynamic Measurements of the Constants of Rochelle Salt,” Phys. Rev., April 15, 
1939, Appendix I. 

’ This method is discussed in a recent paper by W. L. Bond, ‘““The Mathematics of 
The Physical Properties of Crystals,” B. S. T. J., Jan. 1943. 

* The tensor method of writing the elastic and piezoelectric relations is discussed by 
Atanasoff and Hart and by Lawson. See references (2) and (3). 
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transformation of axes, the components of the rotated tensors are eas 
calculated and the elastic and piezoelectric constants for rotated crysta\s 
determined. 

The variation of Young’s modulus as a function of orientation was firs 
worked out by Voigt. In terms of the 7RE angles specifying the orientat 
of a crystal plate, the s compliance modulus (inverse of Young’s Modulus 
is given by the equation 


E E 2 2 me: 2 E - 2 2 
Si1 = Sis(cos Ocos P+ sin’ WY) + (2513 + 544) sin” A cos’ p 


2 2 e 2 - 4 y . . 
X (cos @cos y+ sin” PW) + 533 sin 6 cos‘ p - 2st sin @sinW cosy (A538 


. . 2 . . 
X [3(cos ¢ cos 6cosy — sing sin) — (sing cos 6 cos Wy — cos¢ sin y) 


As discussed in Chapter II by W. L. Bond,” the JRE angles are meas 


ured as follows: Taking the X’ axis along the length of the crystal, the 
along the width, and the Z’ along the thickness, the angle @ is the angi 
between the Z or optic axis and Z’. ¢ is the angle between the projectio: 
of the Z’ axis on the XY plane and the X axis, while y the skew angle i: 
the angle between the length and the tangent to the great circle which co: 
tains the Z and Z’ axes and the length of the crystal XY’. A crystal having 
its thickness along the X axis (X-cut crystal) will have the angles 

6 = 90°; ¢ = 0°; y variable but equal to 90° when the length coincides 
with the Y axis. Under these conditions 


id y + 4 y ° 2 2 
s* = sti sin V+ (2si3 + sha) sin’ y cos 
+ $33 cos* y- 2st, sin’ y cos p 


This equation has been made use of in evaluating the elastic constants 0/ 
quartz as shown by equations (A.20). For this equation Az was measure! 
from the Y axis rather than from the Z as in the 7RE angle and 


A, = yp — 90° (4.60 


Since from equation (A.19) the frequency of a long thin crystal in longitudi- 


nal motion is known to be 
1 7 
f = 4/ ; (A.19 
2¢ pSii 


the longitudinal frequency of any oriented crystal can be calculated from 
equations (A.58) and (A.19). 

It is the purpose of this section to show also that the temperature coetii 
cient of the longitudinal frequency of any oriented crystal can be calculated 
provided we know the temperature coefficient of the six elastic constants 0! 


1 Methods for Specifying Quartz Crystal Orientation and their Determination | 
Optical Means,” this issue of the B. S. T. J. 
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quartz, and that regions of low temperature coefficient crystals can be located 
for this and other simple modes of motion for which the frequency can be 
calculated in terms of the elastic constants. 

Differentiating equation (A.19) with respect to ¢ the temperature 


df_ 1 1 
dt 2 


df 
dt 
f 


where 7, the temperature coefficient of the quantity a is defined as the rate of 
change of a with temperature divided by the value of a. The temperature 
coefficient of the length f = X’ is 7.8 parts per million per degree centigrade 
along the optic axis, and 14.3 parts per million perpendicular to it. For a 
general orientation, the temperature coefficient of length varies as 


E 
dp ds}, 


= Ty = =Ty — UT e+ TH 


| 
J 


Ty = 14.3 — 6.5(sin® @ cos’ y (A.62) 
Since the total mass remains the same when the crystal expands, the tem- 
perature coefficient of the density is the negative of the sum of the coefficients 
of the three axes or 
T, = —36.4 (A.63) 
Hence the temperature coefficient of frequency becomes 


RE’ 
ds}; 
$ a 2 1{ dt 
T; = 3.9+ 6.5 sin 6 cos y — . =i (4.64) 
- Sit 
Differentiating equation (A.58) we have as the temperature coefiicient of a 
general orientation 


T; = 3.9 + 6.5 sin® 6 cos? y 


Si TF, (cos’ 6 cos’ y + sin’ py)’ 4 (2si3T,, + siaT oF, ) xX 
sin’ @ cos” ¥(cos” @cos y + sin’ py) 4 S331 45, 
sin‘ 6 cos* y — 2sts T,£, sin @ sin y cos py X 
[3(cos¢ cos 8 cosy — sing sin)” — (sing cos @ cosy + cos¢ siny Py (65) 
wee Ley CUS | (A.65) 


)| 8 2 Pea, © 2 p\2 E - 8 2 
‘sii(cos @cos y + sin ¥) + (2si3 + say) sin A cos y X 
2 2 mene - 4 4 
(cos @cos » + sin W) + 533 sin @ cos 


| —2sf, sin 6 sin y cos ¥[3(cos ¢ cos 8 cosy — sing sin )° 


— (sing cos 6 cos + cos¢ sin y)’| 
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Hence since the elastic constants are definitely known, the temperatur 
coefficient of any longitudinally vibrating crystal can be obtained wher 
the separate temperature coefficients are evaluated. 

The temperature coefficients appearing in equation (A.65) can all be evalu. 
ated from the temperature coefficient angle curves for X cut rotated crysta\s 
shown by Fig. 1.19. For an X cut crystal equation (A.65) reduces to 


T; = 3.9+ 6.5 cos 
Em 4 eT Em,\.:2 s 
(Sulsf, sin W + (25i37s,, + S4a75#,) sin’ Y cos p 
P . 4 Bon “9 | 
an + 533T,,, cos W — 25147,%, sin’ Y cos ¥ 
lw 2 / E s_% 2 
2 S11 Sin y + (2513 a S44) Sin y cos y | 


+ 533 cos’ y — 2st; sin’ y cos y 


(4.66 


The value of 7,4, is obtained directly for 42 = 0 or Y = 90°, for T; 
and hence 


Tf, = 11.8 (A.67 
Taking three other angles and solving for the remaining constants we find 
T.8,si4 = —5310; — (2s137s,, + StaTs#,) = 45,130; 


(4.68) 


T 535533 = 17,400. 


Inserting the values found for the elastic constants, two temperature coefi- 


cients are determined, and one relation is given between the others, 


TF, = +119; T.,; = 182; TF, — 1112 7,8, = 228.2 (A.69) 


833 844 813 


The valyes of (A.68) are sufficient to determine the temperature coefficient 
of long thin crystals cut at any angle, for inserting these values in (A.65) the 
temperature coefficient for any oriented crystal in longitudinal vibration 


is given by 
T; = 3.9 + 6.5 sin’ 6 cos’ y 
/ +755 (cos? 6 cos? y + sin? y)? + 22,565 sin? @ cos* y 
| (cos? 6 cos? y + sin? Y) + 8700 sin‘ 6 cos4 | 
| +5310 sin 6 sin y cos y[3(cos ¢ cos 8 cos Wy — sing sin py)” 
— (sing cos 6 cosy + cos¢ sin w)?} 
127.9 (cos? @ cos? y + sin? y)* + 175.8 sin? @ cos? p 
| (cos? @ cos? y + sin? y) + 95.6 sin‘ @ cos* y | 
| +89.2 sin 6 sin ¥ cos ¥[3(cos ¢ cos 6 cos YW — sing sin py)? 


| 


— (sing cos 6 cosy + cos¢ sin p)?| 


The only regions of low temperature coefficients are the regions for which 


the two big middle terms are small which requires that @— 0, or y — 0) 
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The first region would be a Z-cut crystal with its length somewhere in the 
YY plane and would result in a temperature coefficient of two parts per 
million negative. Such a crystal] is not of much interest since there is no 
piezoelectric constant for driving it. The other region y — 90° would also 
result in the length being near the XV crystallographic plane, but would 
allow the major surface to be made perpendicular to the X axis and hence 
would allow the crystal to be driven piezoelectrically. By allowing ¥ to 
be slightly greater than 90°, the fourth term in the numerator can be made 
slightly negative and of a value greater than the two positive terms. This 
results in the +5° X-cut crystal having nearly a zero coefficient and this 
angle is the most favorable one for a low coefficient longitudinal mode of 
motion. All other directions have a negative temperature coefficient. 

The remaining temperature coefficients of the six elastic constants can be 
evaluated from Fig. 1.12, and equation (A.22). The frequency temperature 
coefficient can be expressed by the equation: 


T, =3.9+ 65 cos’ 6 


Em - 2 , EB mp 2 on EB. 7 
“ 1| cop 7.8, sin 6 + c437 2, cos 6 + T-.8, ci4 sin 20 A.71) 
E a gE 2 \ E . 
2 C66 SIN O + C44 COS O + Ci4 Sin 26 
since in terms of the 7RE angles the series of crystals is given by ¢ = —90; 


4= 90 — Ao ;~ = 90°. Taking the AT, BT, and Y-cut, whose coefficients 


have accurately been determined, we have 


TABLE VI 
RE 
Crystal Value of Ae Value of 6 Ty 7 of" 
Cut © 6 66 
AT 4+35.25° 54.75 0 —12.0 29.39 X 10!° dynes/cm? 
y 0 90 +86 164.2 40.50 
BT — 49 139 or —41 0 —15.2 O68 86 


From these data and equation (A.71), the three temperature coefficients can 


be evaluated as 


TF, = 164.2; T.#, = 165.7; T#, = +90.2 (A.72) 


To convert these into compliance temperature coefficients we have to 
make use of the relations of equations (A.8) 


Ch zt 

gE E E 44 . E —C14 : 

Seg = 2(si1 — Si2) = : 3 g2? i) Mie a ie gi. ? 
C44C66 — Ci4 2 (C44 Coe — C14) 


- 

7 i C66 

Sua = “38 E2° 
C44 Cog — C1 
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Differentiating these with respect to t, we have 


re Poy 2160. 
ans — 44 C66 Tp ah 14 , 
Vg, = Tok, — | EE Rk? (Tek, + Teh, | + [eee Tf, 


C44. C66 — C44Ce6 — Ci4 


(A.73 
Par cf, 28? 
44.06 TT. Tm 2 “C14 2 
E HE _E2 (7 cf, a 1 ef.) a E E E2 1 ao 


C44 Cog — Cis C44C66 — Ci4 


Py ct 7B? 
~_ 14 C66 nn ane “C14 
E —_ y 4 
1 EE HE _E? (TF, + ra) + EF SS R? TF, 
C44 C66 — C14 C44 C66 — C14 
TABLE VII 


Temperature Present Previous 


Coefficient Determination Determination Bechmann 


T= 11, +12 +11.5 


Sit 

ie . —1265 —1125 
T ~294.8 — 238 —148 
S13 

it +123 +113 
id +213 +180 
533 
Te 
S44 
T= —133.5 —-119 
566 


+189 +175 


TABLE VIII 


Temperature Present Previous Atanasoff 


Coefficient Determination Determination & Hart Bechmann 


Te —46.5 —49.7 —48 


Ci 
: — 3300 35 —3000 | —2115 
1 
T —697 — 5380 —530 
C13 | 
Te +90, +107 482 +110 
T —204.: : —213 
C33 
TE 
Cus 
r® +164. +170.1 +144 +199 


C668 


—165.7 ! — 169 —151 — 199 


Inserting the numerical values for the elastic constants and the tempera ure 


coefficients we have 


T,#, = .8837,%, + .10717,%, = —134.5; 
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The value of 7,#, provides a check on the accuracy of measurement since 
it has been measured in two ways. The agreement is within about 2 per 
cent which shows the probable accuracy of the measurement. Combining 
the coefficients of (A.69) with those given by equation (A.74), the complete 
temperature coefficients are given in Table VII together with previous 


‘ ° 11,12 rn ; ; ~ E . 
determinations '~. The present determination differs from a previous 


determination by the writer due to the use of the elastic constants found 
here rather than Voigt’s constants. 

~~ ° A . *y . 

rhe temperature coefficients of the c; ; elastic constants are easily obtained 
from the s;; constants by employing the relations of equation (A.8). These 
result in the temperature coefficient values for the c constants given in 
Table VIII. 

By using the elastic constant data, the temperature coefficient data, and 
the equations of transformation for rotated axes it is possible to calculate 
the frequency and temperature coefficient of any simple mode for any 
orientation. Examples are given for a face shear mode anda thickness 
shear mode in a previous paper ‘‘Low Temperature Coefficient Quartz 
; i dae ; & 
Crystals.” This paper shows contour maps for low temperature coeffi- 
cient crystals of these types. 

1 The first determination of the temperature coefficients of the writer was given in a 
paper “Electrical Wave Filters Employing Quartz Crystals As Elements,” B. S. T. J., 
July 1934, p. 446. A redetermination using better temperature coefficient data was given 
in a paper “Low Temperature Coefficient Quartz Crystals,” B.S. T. J., Jan. 1940. The 
present determination uses the same temperature coefficient data but slightly different 
elastic constants which results in slight changes in the temperature coeflictents. 

2 A partial determination of the coefficients was made by Koga, Rep. Rad. Research, 
Japan 6, 1 (1934). Other complete determinations are R. Bechmann, Hoch: tech. U. Elek 
Akus. 44,145 (1934) and Atanasoff and Hart, Phys. Rev., Vol. 59, No. 1, Jan. 1, 1941, 
pp. 85, 96. 

B.S. T. J., Vol. XIX, 74 (Jan. 1940). 





CHAPTER II 


Methods For Specifying Quartz Crystal Orientation and 
Their Determination by Optical Means 


By W. L. BOND 


2.1 QUARTZ AND ITS AXES 

The chemist describes quartz as silicon dioxide, SiOQ:, crystallized jy 
hard, brittle, glass-like, six sided prisms, often with pyramidal terminations: 
melting point 1750° Centigrade, density 2.65, hardness on Moh’s scale 7. [| 
transforms from alpha to beta quartz at 573°C under atmospheric pressure. 
Under stress it transforms at lower temperatures. Alpha quartz is in- 
soluble in ordinary acids but soluble in hydrofluoric acid; and in hot alkalis, 

At first glance we might say that it had hexagonal symmetry but if we 
etch two adjacent pyramid faces we, find that the microscopic etch pits 
are of different shape, hence the faces cannot be equivalent. It has three 
axes of two-fold symmetry and one axis of three-fold symmetry. Let us 
also remark that it does not have a center of symmetry or a six-fold axis 
Figure 2.1 shows us that the three two-fold axes are perpendicular to the 
three-fold axis and are 120° apart. If they were not like this, they would 
not be self-consistent. 

As we examine more and more quartz crystals we find that there is a 
tendency for pyramid faces to be alternately large and small, the larger 
faces being brighter than the smaller faces. Also the etch pits of alternate 
faces are similar. (The etch pit study is a powerful tool in determining 
crystal symmetry.) Further, two other “kinds” of faces are quite 
commonly found. If we draw such a crystal as though equivalent faces 
were of equal size we get such a picture as Fig. 2.2. It is an idealized figure 
used to illustrate the symmetry of quartz. The prism faces are marked 
m, the six faces marked r ‘“‘constitute the primary rhombohedron”’—the 
ones we called the large bright pyramid faces. The crystallographer 
thinks of these six faces as pieces of the faces of a rhombohedron. (A 
crystallographer’s rhombohedron is like a cube stood on one corner, then 
the opposite corner pushed in a little towards the other, or pulled away 
from it. He thinks of it always as standing on this corner, as Fig. 2.3. 
The z faces constitute a second rhombohedron—the secondary rhom- 
bohedron or minor pyramid faces. The s and x faces illustrate a further 
property of quartz. Figure 2.3 differs from its mirror image so that we 

224 


migh' 
one’s 
right: 
in Fis 
has § 
these 
its oc 

Th 


one k 


pare 
in t 
and 
ther 
posi 
han 

L 
side 
the 
Thi 


a Pp 





SPECIFYING QUARTZ CRYSTAL ORIENTATION 


might expect to find two kinds of quartz that are related to each other as 
one’s right hand is related to his left. We do find them and call them 
right-hand quartz and left-hand quartz respectively. They are illustrated 
in Fig. 2.4. These conventional] figures are shown in many texts but no one 


has seen such perfect quartz crystals. They are drawn possessing just 


these faces and no others merely to illustrate the symmetry of quartz and 
its occurrence in right-handed and left-handed forms. 

These figures are also useful in defining how a blank shall be cut from 
one kind of quartz. It is found that if a crystal be compressed with forces 


A3 








A3 


Fig. 2.1—Hexagonal axis system 


parallel to a pair of sides of the hexagon an electric polarization takes place 
in the direction of the forces. The edge “modified” by the presence of s 
and x faces becomes negative. If we allow these charges to leak off and 
then suddenly release the mechanical forces the “‘modified’”’ edge becomes 
positively charged as the crystal expands. This is true for both right- 
hand and left-hand crystals. 

Let us now talk about right-hand quartz. Since expansion is con- 
sidered as a positive strain (contraction as negative) it is decided to take 
the positive end of the electric axis as pointing towards the modified edge. 
This gives us a positive charge at the positive end of the electric axis when 
a positive stress (tension) is applied along this axis. This positive direc- 
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tion of an electric axis is taken as the positive x axis of a right-hand 
The z axis is taken along the axis of th; 


rectangular coordinate system. 
hexagonal prism, and since the x axis is an axis of two-fold symmetry 














Fig. 2.2—An idealized quartz crystal 


can take either end of the prism as the direction of +z. We then choose ) 
to form a right-handed coordinate system. (In a right-handed system if a 
right-handed screw turns about the z axis in the sense x to y it would ad- 
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vance in the positive z direction.) The y axis will always lie directly under 
a major rhombohedral face. 

We could define the x, y and z axes for a left-hand crystal as forming a 
left-hand system. Though this is a useful conception in mathematical 


A rhombohedron 


LEFT HAND RIGHT HAND 
QUARTZ QUARTZ 


Fig. 2.4 


studies, we can dodge this “‘double standard” by a simple device. For use 
as a crystal circuit element, left-hand quartz can be used just as well as 
right-hand quartz. In designing such an element it suffices to think al- 
ways in terms of right-hand quartz and issue specifications for this kind 
only, using always right-hand coordinate systems. If now for left-hand 
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Simple crystal cuts 
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Fig. 2.6—Simply rotated cuts 
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crystals we mark the negative end of the electric axis as positive we can 


treat it exactly as though it were a right-handed crystal. 

The first plates used were x and y piezoids (squeezing solids). For 
these simply described cuts one does not need to know the quartz “ handed- 
ness.” These crystals had large frequency temperature coefficients. 
But when Lack, Willard and Fair brought out the low temperature coeffi- 
cient AT plate, its more complicated orientation required the right-left 
differentiation. The AT, the subsequent BT, CT, DT, etc., were thought 
of as y-cuts rotated through various angles about the edge that coincided 

i 


with x. For example, the AT was a +35{° cut, or was a y plate rotated 


353° about x; the BT was a —49° cut. Their orientations are illustrated 
in Fig. 2.6. 

As more complicated orientations were designed to give even better 
temperature coefficients at extreme frequencies the description became more 
difficult, requiring the specification of two or three angles. Many schemes 
for specification have been devised but the Institute of Radio Engineers is 
recommending the adoption of a system we shall call the I.R.E. system. 

The crystal designer has the problem: “How shall I orient the length, 
width and thickness of a piezoid with respect to the x, y and s axes so as to 
give the required electrical properties? He thinks in terms of fixed x, y 
and z axes, variable piezoid edge directions. The crystal cutter has the 
converse problem. “How shall I orient the x, y and z axes of the crystal 
so that fixed saws will give the required surfaces?”’ For this reason the 
most convenient orientation angles from the designer’s viewpoint may 
not be the simplest from the cutter’s viewpoint. Also the translation 
from one set to the other may not be simple. 

The early methods of orientation specification were somewhat chaotic. 
There was no overall plan of what angles were to be specified and from 
what axes they were to be measured. Each group of crystals was a law 
unto itself. 


Tue I.R.E. ORIENTATION ANGLES 


The relations between the x, y and z axes of the crystal and the length, 
width and thickness of the piezoid are given in Fig. 2.7. 

The position of Fig. 2.7 may be considered as a result of turning the piezoid 
through the successive angles ¢, 6, y starting from an initial position length 
parallel to x, width parallel to y and thickness parallei to z as in Fig. 2.8. 
First the crystal is turned through angle ¢ about z in the direction shown 
in Fig. 2.7. Then it is lowered through angle @ about an axis parallel to the 
width direction, again in the direction shown in Fig. 2.7. Finally it is skewed 
through an angle y about an axis parallel with the thickness direction in 
the sense shown in Fig. 2.7. 
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Fig. 2.7—-The I.R.E. orientation angles 
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Fig. 2.8—The initial position 0, 0, 0 for the I.R.E. ang.es 
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Tue I.R.E. ANGLES FoR A FEW STANDARD PIEZOIDS 


Name o 9 J 
Z 0 0 0 
x 0 90 90 
— 18° Filter 0 90 108 
+5° “ 0 90 85 
y 90 90 90 
AT —90 543 90 
BT —90 —41 90 
ay —90 52 90 
GT —90 38° 52’ +45° 
MT 6° 40’ 50° 28’ 79° 36’ 
NT a 40° 40’ 77° 40’ 


2.2 ORIENTATION BY NATURAL FACES 


With well faced material one can do an accurate job of orienting without 
X-rays if he knows the faces of quartz thoroughly. 

The quartz rhombohedral faces are highly perfect and polished, the 
major often more so than the minor. With two such faces a device like 
that illustrated in Fig. 2.9 could be used to give an orientation accurate to a 
minute or two. An adjustable base, symbolized here as a ball and socket, 
is adjusted so that the eye centers the lamp filament image on the cross 
hairs, first for one face, and then, turning the base about on the reference 
table it is adjusted for another face. When the images all pass through 
center as the base is turned on the table the optic axis is perpendicular to 
the table. When any one image is centered, the electric axis is perpendicu- 
lar to the plane of the paper. This with the already mentioned fact that 
AT plates are cut near a minor face and BT’s near a major, allows us to 
cut the crystal accurately. 

Although the rhombohedral faces are highly perfect the prism faces never 
are. On the prism face, true prism faces that are very short in the z direc- 
tion alternate with short rhombohedral faces to give the general contour 
aslant. These ‘‘steps’’ give the face a striped look. The stripes are known 
as growth lines or striations. They are parallel to x and can be used to 
find x to a degree or so. If we sight on striations on two adjacent faces 
we can locate the optic axis to nearly the same accuracy since the optic 
axis is perpendicular to the striations on all faces. 

There are several indications that help us find, from the prism, where 
the major rhombohedron would be in the absence of such faces. Some 
crystals grow in the form shown in Fig. 2.10. They are symmetrically 
doubly terminated and a very narrow prism is found under the major 
thombohedron, a wide face under the minor. Hence given a portion like 
that enclosed in the dotted line we could deduce the complete orientation. 

Some crystals grew out at right angles to a wall and because they grew 
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Fig. 2.9—Optical orientation by reflection of light from faces 











Fig. 2.10—A type of quartz growth 
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along a 3 axis in one direction only, the x axis does not appear to be one of 
two-fold symmetry. Such a crystal is illustrated in Fig. 2.11. Here the 
prism faces under a major rhombohedron are tapered and bright, the 
prism faces under a minor are relatively parallel sided and very dull. The 
bright prism faces are much more nearly parallel to the optic axis than the 
dull ones. Again, given a portion of the prism we can deduce the orien- 
tation. 





Fig. 2.11—Another type of quartz growth 


i) 


3 Fiat Lay Cuttinc 

Flat lay cutting takes advantage of the fact that, although tapering 
quartz prisms have their faces non-parallel to the z axis the prism faces are 
parallel to the x axis. A crystal is cemented prism face down, to a mount- 
ing plate. The mounting plate is tilted and turned on a base plate to render 
the optic axis parallel to the long edge of the base plate. This is done in 
a conoscope. Now the edges of the base plate are the x, y, 2 axes of the 
crystal. ; 

The crystal can now be cut directly into wafers for dicing into AT’s, 
BT’s etc. by mounting on an angle bracket as shown in Fig. 2.13 or cut into 
X sections from which AT or BT bars shall be made by merely sawing down 
the length. Again z sections can be cut by cross cutting. Good z sections 
can be made in this way in the total absence of faces. These sections can 
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Fig. 2.13—Direct sawing of the slabs - 
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then have their x axes determined by etching and X-rays and cut up by the 
Z section cutting method. 
By turning the base plate on the angle bracket and dicing the wafers at 
an angle any orientation can be obtained. 
2.4 Z SECTION OR VERTICAL CUTTING 
Javing a true Z basal section it is first marked for the +x axis. The 
simplest procedure is to use the star mark; for right-hand quartz (R.H.Q. 


LHQ RHQ 
+ + 


Fig. 2.14—Marking the ‘“‘sense’’ of righthand and lefthand sections 
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Fig. 2.15—Making the rotations A: 


the rays should point toward the plus electric axis, for left-hand quartz 
the rays should point towards the negative electric axis. 

The section is now placed on the carriage plate, one ray pointing towards 
the operator (which ray is decided on the basis of the economy of quartz). 









































Fig. 2.16—The rotation A» 
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Fig. 2.17—The slab after the rotations Ai and A, 
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[he section is then rotated clockwise on its base, through angle A; as in 
Fig. 2.15 and cemented in this position. 

The carriage plate is then transferred to a diamond saw angle bracket of 
tilt de, as in Fig. 2.16, and the crystal is sawed into slices slightly thicker 
than the required final thickness f. 

The operator turns these slices down flat on the table of a dicing saw as 
in Fig. 2.17 by rotating the slices 90° clockwise about the axis AA, then turns 

F the slice through angle A; as in Fig. 18 and makes a cut. The plate is fin- 
ished as shown in Fig. 2.18. 

Since the angle bracket is not reversible, negative A» angles are cut by 

adding +180° to A; and reversing the sign of A3. 


Tue A ANGLES FOR SOME STANDARD PLATES 


Cut A, Ao Ai 
x 90° 0 0 
y 0 0 0 
z 0 or 90° 90° 0 or 90° 
—18 90° 0 + 18° 
oA§ 90° 0 —5 
AT 0 353 0 
BT 180° +-49° 0 
Ci 0 38 0 
DT 180° §2° 0 
GT 0 Si° 7° +45° 
MT 96° 40’ 39° 32’ —10° 24’ 
NT 99° 25’ 49° 20’ — 12° 20’ 


2.5 THE RELATION BETWEEN THE I.R.E. ANGLES ®@) AND THE Z SECTION 
ANGLES A,, Ao, Az 


It can be shown that: 


Ay = 90 +o 
Ag = 90 7 6 
Aj; = — 90+ y 


2.6 POLARIZED LIGHT AS APPLIED TO CRYSTALS 


Light consists of electromagnetic ‘‘vibrations.”” The vibrations are per- 
pendicular to the direction of propagation but ordinarily helter-skelter in 
all directions perpendicular to the propagation. The color of the light is 
determined by the vibration frequency, blue vibrating more rapidly than 
red. Ina vacuum, light travels at 186,000 miles per second (3 X 10'° cms 
per second) all colors at the same velocity. On entering a transparent me- 
dium the velocity is reduced, ordinarily blue being slowed more than red. 
The frequencies are unaltered on entering the medium. 

Light traveling through a uniaxial crystal in the direction of Fig. 2.19 
breaks up into two components that travel at different velocities. For one 





238 BELL SYSTEM TECHNICAL JOURNAL 


of these components the vibration is all in the plane of poz, of the 
the vibration is all perpendicular to the plane of poz. 
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Fig. 2.19—The velocities of light in a uniaxial crystal 


A plot of the propagation velocities for all directions is a surface of two 
sheets, one octant of which is shown in Fig. 2.19. One sheet is a sphere; 





f two 
here; 
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the other sheet, which is an ovaloid of revolution, touches the sphere at the 
two points where the double sheet is pierced by the optic axis. Of the two 
ravs traveling along one line, the one with a velocity corresponding to the 


sphere 1s called the ordinary ray; the one with a velocity corresponding 


to the ovaloid is called the extraordinary ray. For quartz the ovaloid is 
prolate and lies inside the sphere. For tourmaline the ovaloid is oblate and 
ies outside the sphere. The small arrows show the direction of vibration 
Each of the components is said to be polarized since for each all the vibra- 
tion is in one direction. 

Since both sheets are surfaces of revolution with the optic axis as the 
ixis of revolution, we can never tell the x axis from the y axis by optical 
means.* Only the s axis is a unique direction and can be determined op- 
tically. If this figure is taken to represent the case for blue light there will 
be a slightly larger but similar figure for red light since, in the crystal, red 
light travels faster than blue light. 

Some kinds of crystals have velocity plots for which the double sheet 
surfaces touch at four points. Hence they have two optic axes and are 
called biaxial. All hexagonal, rhombohedral and tetragonal crystals are 
uniaxial, all others except the isometric ones are biaxial. Rochelle Salt is 


biaxial. 
2.7 POLARIZERS AND ANALYZERS 


In the Nicol prism means were found to eliminate the ordinary ray; the 
other is transmitted as polarized light. That is, ordinary light of any or 
all colors upon passing through a Nicol prism emerges as plane polarized 
light with no change in color. 

Transparent colored media appear colored because they absorb some 
colors of light more than other colors. In cplored crystals the two rays 
themselves often differ in their color absorption so that the crystal as viewed 
by means of the ordinary ray seems to be of a different color than as viewed 
by the extraordinary ray. Quinine iodo-disulfate, or Herapathite, absorbs 
most visible colors of one ray almost completely; transmits about 60% of the 
visible colors in the other ray. Hence light emerging from this crystal is 
almost completely polarized. In the commercial product called “polaroid”, 
myriads of such crystals, microscopic in size, are contained in a celluloid- 
like sheet and oriented by stretching the sheet. This material now re- 
places Nicol prisms for all but the most exacting uses. 

If we put two identical polaroid sheets together with their transmission 
vibration directions parallel as in Fig. 2.20 we can see through them but if 
their transmission vibration directions are at right angles we cannot see 

* Methods depending on etch pits are excluded. They are optical only in the sens¢ 
that observing natural faces is optical. 
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through them because the second sheet can transmit none of the vi! 
transmitted by the first sheet. As we rotate the second sheet bac fro; 
complete extinction to “best transmission” the transmitted light increase 
sinusoidally. In any such arrangement as Fig. 2.20 the first sheet is cajjp 
the polarizer, the second is called the analyzer. The name analyzer \ 
chosen because light that can be extinguished by means of a suital)ly ; 
tated analyzer must be plane polarized, and it must be vibrating at righ: 
angles to the transmission vibration direction of the analyzer when set jo; 
extinction. The transmission vibration direction of a polaroid plat 
hereafter be called its vibration axis. 

Let us go back to Fig. 2.19 and cut out from around the point f, the sma 
tangential crystal plate shown magnified in Fig.2.21. Here p is the direct; 
of propagation as before, and z is parallel to the optic axis. Also s whi 








UNPOLARIZED y} / POLARIZED 

OR PLAIN LIGHT VIBRATING 
LIGHT VERTICALLY 
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Fig. 2.20—Light polarization 














in the plane of p and g, is the direction of slow vibration while f whic! 
perpendicular to this plane is the direction of fast vibration. The vibra 
~ tion frequency is really the same for both. ‘Slow vibration’? means ‘\ 
bration direction for slow transmission.” All directions of propagat 
that have this vibration axis have the same velocity. 

In Fig. 2.22 we have placed this plate between “crossed polaroids’’— tha’ 
is polaroids set for extinction. The slow direction makes an angle a wi! 
the polarizer vibration axis. When this vertical-polarized ray of intensit 


/ enters the crystal it breaks up into components, one of intensity: 
I cos a vibrates along s 

and one of intensity 
I sin @ vibrates along f 


as illustrated in Fig. 2.23. 
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2.21—Light breaks into components in the crystal 


POLARIZER ~ ANALYZER 


Fig. 2.22—Transmission when a crystal is placed askew between crossed polarizers 


If a = 0 the fast component reduces to zero and the slow component goes 


through the crystal unchanged hence emerging plane polarized. It can 
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then be extinguished by the analyzer. Ii a 90° the slow comp 
reduces to zero and the fast one goes through und hanged and again ( 
extinguished by the analy zer. This effect can be used to check c1 
orientations. Such an instrument fitted with a divided circle used to 1 
ure @ is called a stauros« ope. The stauros ope often uses a sper ial ana 
capable of better determination of extinction setting 

If @ is not 0 or 90° two components traverse the crystal and recomb 
the boundary. These two components are at right angles to each « 
they are of unequal intensities, and they differ in phase because they 


eled at different speeds. 


Tp 





Ps 










ponents from | ig 22 


Now v, and v, have the same frequency / so that in unit time each mak 








F wave-lengths. This requires that the slow ray have F wave-lengths 


distance v, and hence that each wave havea length: 
Similarly 


. - l 
In a distance / there are | fast waves and " slow ones. Let us say 


‘Af / 


there are V more fast waves than slow ones in the distance ¢. Consequent 


: t t ; . : oP P 
N = ~- which, from (2.3) and (2.3’) may be written: 
i & 
ae | tl 
N= . } 
Uy t 


All the data on light are given in terms of wave-lengths in a vacuum, 


in terms of frequency, so we will assume that in a vacuum the wave-lengt 
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this light is A, and as in a vacuum its velocity is | 


second) an equation similar to (2.3) would tell us that A j and hence that: 
y 

I (2.5) 
A 


With (2.5) we can rewrite (2.4) as 


Phe ratio of the velocity ina vacuum to the velocity in a medium is called 
the refractive index of the medium commonly given the symbol ». For 
most transparent materials 2 is between 1.3 and 1.8 


\e write these refractive indices as 


V V 
= my and = Mm, respectively. 
vs Us 
Now (2.4) becomes: 
a se 
N : . (ny — Nz) (2.7) 


After passing through the crystal plate of thickness /, Fig. 2.21, the two 
ight components recombine. They are polarized at right angles each to 
each; they are of unequal intensities, and they differ in phase by N wave- 
engths as given by equation (2.7 

If the crystal were vanishingly thin the two components that recombine 
would be effectively in step or in phase. In Fig. 2.24 we have plotted vibra- 


tion amplitude against time for these two components. They are sep- 


arated for clarity. In the upper curve the slow vibration is shown as ver- 
tical, in the center curve the fast one is shown as horizontal. In the lower 
curve corresponding points have been added vectorially. From actual con- 


struction we see that in the resultant curve the vibration is always parallel 
to the line 44’. Hence it is plane polarized and can be extinguished by 
means of an analyzer. 

Let us now consider the case of a thicker plate for which the slow wave 
gets a quarter-wave-length behind the faster one. This case is plotted in 
ig. 2.25 in the same way that the previous case was plotted in Fig. 2.24 
It turns out to be a space curve like a slightly flattened corkscrew. Viewed 
ilong the axis it looks like an ellipse, as shown to the right of the space 
curve. If the slow ray had lost but a little with respect to the fast ray, we 
would have gotten a very flat ellipse. If the two components had had the 
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same amplitude with the quarter wave phase difference the end view in fz 
2.25 would have been a true circle. Now since the vibration, in these casvs 
is not all in one plane, the light cannot be extinguished by an analyzer, i: 
is not plane polarized light. In the one case it is called “‘elliptically po 
ized”’ light and in the other, ‘circularly polarized” light. If the slow 
loses an integral number of wave-lengths, it makes no difference; only fr 





Fig. 2.24—The recombination of the light components after passing through a 
thin crystal 


tions of wave-lengths count, except that if several wave-lengths are present 
a thickness that is right to give an integral number of wave-lengths for one 
color may give an integer plus a fraction for some other color. If the thick- 
ness is fairly small, this may cause spectral colors from white light. For 
thick plates the wave-lengths so overlap that the field appears colorless 
but dark or white according to the value of the angle a; i.e., if @ is zero 01 
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40° the field is dark; if a + 0 or 90° the field is bright. In Fig. 2.26 we illus- 








in Fig 

© Cases trate how, for one color of light the polarization starts out as plane at the 
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thick- Fig. 2.25—Recombination more generally 
For 
lorless crystal boundary, passes through elliptical to circular polarization, then flat- 


ero or tens out the other way through elliptical to plane polarization at a distance 
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in the crystal corresponding to the slow wave being one wave-length be 
the fast one. 

In Fig. 2.27 we show two AT plates resting on a reference surface. | 
first crystal the propagation is perpendicular to x, in the second crystal! 
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Fig. 2.26—-How the kind of polarization changes with crystal thickness 
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Fig. 2.27—The normascope principle 


along x. If the reference surface is the reference table of a simple stau: 
scope the edge of the first crystal will appear dark because a = 0° or 0%) ; 
the edge of the second crystal will appear bright because a@ is not 0° or 
Actually this a for an AT plate can be +35° + any multiple of 90° becaus 
we don’t know whether z stands out to the right or to the left. Hence 1 
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1 be 


reading might be, for instance, —35°, +35°, +55°, +125°, et This is 
the principle of the normascope used to identify the x direction for crystal 
; adjustment. 
Stal 


Let us study these relative phase shifts at different angles near the optic 
axis. Now quartz has an optical complication beyond that just described 
it rotates the plane of polarization of plane polarized light traveling along 
he optic axis. This complicates our present attempt to build up a back- 
ground sufficient for an understanding of the conoscope. But the cono- 
scope finds the optic axis for other crystals that do not rotate the plane of 
polarization, tourmaline for example; so we will ignore this rotation, to be- 
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Fig. 2.28—How the polarization changes with propagation direction in a crystal plate 


gin with, in order to arrive quickly at some useful conclusions. We will 
later explain how optical rotation modifies these conclusions. 


Consider then the crystal z section shown in Fig. 2.28. A source s sends 


e monochromatic light through the polarizer which passes only vertical 
vibrations. We will assume that the light passes in and out of the crystal 
without a deviation of path. Since the vibration is in the plane of 2 and 
p; (its direction of propagation) the ray does not break up inside the crystal 

Sstaul but is propagated as plane polarized light, unchanged. An analyzer set 

or ) for vertical extinction could then extinguish this ray. 

or This is true for propagation from s anywhere in this vertical plane 

atoeialt \lso since the vibration is perpendicular to the plane of z and py» the ray 


nee th »» does not break up inside the crystal but passes through and out un- 
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changed. All rays from s in the horizontal plane emerge plane polari ed 
and can be extinguished by an analyzer set for vertical extinction. {he 
ray po is in both these planes so it can be similarly extinguished. 

With the ray p; the situation is different. Here the vibration is not in the 
2 ps plane so the ray breaks up inside the crystal into two components which 
travel with different velocities and recombine in or out of phase to give the 
various degrees of elliptical polarization (including plane and circular 
Hence, an eye looking back along fs, through an analyzer set for vertica! 
extinction, will see light or dark depending on the phase shift VN. Now this 
phase shift for a given thickness of plate is zero along po but increases as 
6 increases (without changing ¢; see Fig. 2.28), passing through one integra 


o=0° 








Fig. 2.29—A plot of phase as a function of ¢ and 6 


value after another. Therefore, as we allow @ to increase, the eye should 
see alternate dark and bright regions. Moreover, since the crystal] is op- 
tically symmetric about z, if @ is changed without changing @, the apparent 
brightness will not change (except that if @ = 0, 90°, 180° or 270° the field 
is dark as we previously explained). Consequently, if we could see along 
all directions at once we would see a pattern of concentric dark rings on a 
dark cross as shown in Fig. 2.29. 

But we can see along all these directions at once if we employ a proper! 
placed lens for a lens can converge all these rays to one point where an eye 
can be placed for viewing. 

Thus an eye at e, Fig. 2.30, will see, in the direction e po’, the ray that 
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started along s po. It will see along e p,’ the ray that started along s py’. 
Every point on the lens will have associated with it a different direction 
in the crystal. Therefore the eye will see a pattern like that of Fig. 2.29. 
[his is the principle of the conoscope. In the conoscope (the name means 
“conical viewing”) the source s is replaced by the image of a source, the 
image being cast by a lens; see Fig. 2.31. Thus by the use of two similar 
lenses we get twice as much working space as one lens would give. 
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Fig. 2.30—The principle of the conoscope 
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Fig. 2.31—A practical conoscope 


Figure 2.32 shows a cross-section of the Western Electric conoscope. The 
graduated dial shaft goes out through the bottom of the tank to give more 
working room-—older instruments had the shaft overhead and it was in the 
way. The light source is a mercury arc lamp with filters to isolate the 
5461A line. The lenses have a converging power corresponding to f:0.6. 
The focus is not changed by changes in the refraction of the oil—in fact, the 
focus is the same with no liquid in the tank as when filled with liquid. This 
is of some interest for those who might wish to use the instrument for 
Rochelle salt and accordingly use a fluid of refractive index about 1.495 
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nstead of the 1.546 of quartz. The dial is graduated into degrees and a 
vernler allows readings to tenths of degrees. 
Ii a crystal plate is held against the glass reference surface one may read 


he angle between the optic axis and the surface normal. One should 


ccasionally check the instrument (against slippage of the dial) by reversing 
the crystal and recentering the pattern. Ii the readings are not identical, 
he dial should be adjusted till they are, ven if the readings are not iden 
tical the mean value should be correct. If one is using the method of ring 























€ IN DEGREES 
































centering, the correction for eccentricity from Fig. 2.33 should be applied to 
this mean value. 

The carriage may be slid back and forth and for very small crystals the 
carriage should be placed so that the crystal is near the center of the tank 
so that very little of the light cone by-passes the crystal. By the use of a 
block a thin crystal can be examined by viewing through its edge or length. 

The carriage can be removed and ‘‘raw”’ crystals examined. The optic 
axis is plainly visible and quite accurate orientations can be made if there 


is not too much opaque material on the crystal. Excessive optical twinning 
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makes a confused pattern but good orientations can be made anyw 
“raw’’ crystal can be mounted adjustably in a jig that is lowered i: 
conoscope, the optic axis lined up, the jig transferred to a saw, and s 
sawed directly. 

Let us turn now to the quantitative analysis of the ring pattern s 
the eye piece when examining a uniaxial crystal. We wish to know t 
of the smallest ring in the field, or rather the corresponding angle 
crystal. This first dark ring (analyzer and polarizer crossed) is the : 
of the slow wave falling one wave length behind the fast one. If th 
thickness (Fig. 2.34) is ¢’ the path length in the crystal is 
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Fig. 2.34—The angle of the smallest ring 


This is to be substituted in Eq. 2.7, namely: 
‘ t 
N= . (ny — ns) 


Now it can be shown that, quite accurately, at the angle @ from the op 
axis: 

= .« 9 (7 ¢ 

Ng — ne = .00917 sin’ 6 [2.9 


where .00917 is the difference in the refractive indices for the ordinary r 
and the extraordinary ray for green mercury light traveling at right angles 
to the optic axis. (These are generally given the symbols , and n, or ) 
and , respectively.) 


t call 
N; = X .00917 sin 6 = 1 
* ~ ) cos 8 
and since \ = .000546 mm. for green mercury light this may be writt 


t’ sin 6 tan 6 = 0.0595 mm. (2.10 
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uw 
w 


vhence we solve for the values in this table 


convergence @= 5 iv 2 .3¢ 
thickness =78 1.94 0.48 0.21 


his shows that if we wish to examine thin plates in a conoscope the lenses 
must be strongly convergent. The conoscope used in the Western Electric 


as a convergence corresponding to about the 20° entry of the table so it 
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Fig. 2.35—-Ring eccentricity 


can be used on crystals down to a half millimeter thick—that is for orien- 
tation studies. In determining handedness we remember that this is a ques- 
tion of the rotatory power of quartz causing the rings to expand or contract 
on rotating the analyzer. Also we said that this rotatory power effectively 
disappears at 15° from the axis. If no ring is found within 15° of the axis 
there is no ring capable of expanding or contracting and we cannot test the 
handedness of such a thin crystal no matter how strong a lens we employ. 
We can then fall back on the succession of colors shown when we rotate the 
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analyzer using white light. Red, yellow, blue are observed for clock 
rotation with IRE right-hand quartz. 

If the s section is not a true one, Eq. 2.9 will be replaced by one allowi: 
this error 6: 


n, — ny = 00917 sin? (8-3) 


This will make the rings non-concentric and although the ‘‘cross”’ j 
section is still the true optic axis the ring centers are not. 

Hence if we are tempted to find the optic axis by centering a large s 
ring in preference to a small fuzzy one we find this eccentricity error must |) 
allowed for. 


2.8 ROTATION OF THE PLANE OF POLARIZATION 


If we have a polarizer and an analyzer set for extinction (Fig. 2.36), the 
place a thin s section of quartz between them, the field brightens up 
can be extinguished again by rotating the analyzer, Fig. 2.37. For the: 
eralogist’s right-hand quartz the analyzer must be rotated 21.7° (y 
light assumed) clockwise to re-establish extinction, counter-clockwise 21.7 
for left-hand quartz. The rotation is more for blue light, less for red. i 
the section is not a perfect s section the rotation is less than this, effective 
disappearing at about 15° from the optic axis. 

A thick slab can be examined ir this Way and, due to the color differ 
in rotation, ‘‘rainbows” will be seen in the quartz when held at just 
right orientation. These rainbows will follow the contours of the specim: 
unless both right and left quartz are present in one piece. When th 
the case the one kind generally occurs as spike- or blade-like intrusior 
the other. It will then cause the rainbows to have sharp, jagged outli: 
bearing no relation to the specimen contour. 

Also, since red, yellow, blue, are here in the order of increasing rotatiot 
if we rotate the analyzer clockwise for right-hand quartz (I.R.E. RHQ 
will pass through best transmission for red, best for yellow and best 
blue in that order so that the field will assume these colors in this order 

With uncut stones this examination is best made under an immersion fluid 
The inspectoscope is made for this work. We spoke of the rotation of th 
plane of polarization and its complicating of the issue for the conoscop 
Due to this the field at the center is not dark when the analyzer and polar 
izer are crossed. Also if we rotate the analyzer clockwise the rings of th 
pattern either expand or contract according to whether the crystal is right 
hand quartz or left-hand quartz (IRE definition). 

A different kind of pattern is visible in the conoscope when viewed px 


pendicular to s, a double set of hyperbolae as shown in Fig. 2.38. This pat 
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D¢ S tern has been used to check the orientation but much grief has ensued due 
to not recognizing one of its properties. This property is that, if z does not 
V1 lie parallel to the crystal boundary the center of the pattern is not per 
pendicular to the optic axis and a rather involved correction must be used. 
This correction reduces the actual angle to about half the observed value. 
[his conoscope is an immersion instrument. ‘The tluid is chosen to have 
7 an index of refraction to match the “ordinary” one for quartz. When this 
is done light is not bent in passing between fluid and quartz. When the 
sie fluid does not match there is a bending and all readings are subject to a 
nu 
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ast Fig. 2.36—Crossed polarizers Fig. 2.37—Rotation of the plane polarizer 
er 
 flui correction. For example if we measure the angle of an AT plate in a fluid 
of the that is too low by .0048 (since my for quartz in green mercury light is 1.5462 
scope this fluid has n = 1.5414), we will get a reading that is too high by a quar 
pola ter degree (the 35° angle will then appear low). A temperature rise of 12° C 
of the will lower the relative refractive indices by this amoun 
rig Also the more perfectly the fluid matches, the more nearly will the rough 
quartz surface disappear and seem smooth and clear. This greatly en- 
d pr hances the sharpness of the rings. 


The refractoscope (Fig. 2.39) was designed by G. W. Willard to tell when the 
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match is good. It uses the elimination of the bending as a test for refrac! 
match. It also demonstrates the existence of two velocities in quartz, 
two images are seen where a glass prism would cause but one. Also 
viewing through an analyzer we see that the two images are caused by p! 
polarized light, the polarization planes being mutually perpendicular. 
If the fluid has an index lower than that of the prism the rays will be 
towards the base of the prism. For this reason, light that reaches the ey: 
from s must travel by the path s gopPoe for the ordinary ray, s q.p.e for 


























Fig. 2.38—Conoscope view normal to the optic axis 


extraordinary ray. Hence the slit as seen through the prism will appear at 


I » e 4 . ° . 
so for the ordinary ray, at s, for the extraordinary while the slit as seen 


directly alongside the prism will appear at s. If the fluid index matches 
the quartz prism for the ordinary ray this ray will be unbent at po and qo s 
that so will appear as a continuation of s. 

If the fluid index is too high the image so will appear to the left of s wit! 
s, still to the right of so. 

As the refractive indices of quartz for the ordinary and the extraordinary 
ray differ by .009 the apparent separation of so and s. represents .009 and 
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can be used to judge the difference between the liquid index and the quartz 


dinary index. 

The liquid can be adjusted by placing a cap face of a crystal on the refer- 
ence glass and setting the dial to read 51.8° the liquid being then blended 
to center the pattern. 

If the refractive index of the fluid is low by an amount L the observed 
reading R must be corrected by adding to it an amount e. 

Since 

(N, — L) sin R = N sin (R—-e€ 


we can compute e«. The accompanying refraction correction nomograph 
was computed from the above equation. If we know that the fluid in the 
conoscope is high, by an amount //, and we wish to know the correction to 
be applied to a conoscope reading R we locate H on the diagonal line HL 
and locate R on the horizontal line R. We join these two points with a 
straight line and read the scale e—e where this straight line crosses the 
curved line e—e. This value tells the size of the correction and whether 
to add it to R or subtract it from R. Conversely, if we wish to find how 
closely the index of the fluid must be held in order to have the correction 


less than say }° at a reading of say R = 50° we join the points R = 50° and 
e = + 2° and find H = + 0.005. A ten-inch-long lucite strip with a 


straight line ruled on it is a convenient tool with which to read this nomo- 
graph. 

We now inquire as to whether the refraction correction can be made to 
annul the ring eccentricity correction. In the appendix it is shown that 
this is done if H = —.530 tan? d where 2d is the distance between the ver- 
tical reticule lines. 

Experimentally it is easy to achieve this balance by using a cap face 
parallel slice. With the cap face against the reference glass and the dial 
reading 51.8° the fluid is blended to make a single ring tangent to both reti- 
cule lines. When this is done for d = 10° the fluid should have a refractive 
index of 1.5228 and the residual errors should be less than 2 minutes for R 


not over 60°. 
2.9 IMMERSION FLUIDS 


In order to match the refractive index of quartz we blend a substance 
which has an index that is too large with one that has an index that is too 
small. Such blended substances should be liquid at room temperature and 
hence should be perfectly mutually soluble. They should have low vapor 
pressure so that they do not evaporate quickly and should be harmless to 
the operator. Also they should be nearly colorless and clear. They should 
be fluent enough to be easily drained from the crystal and should have a 
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sufficiently high flash point that they would not present a fire hazard. 

odor should not cause distress and finally the cost must be reasonable. 
Dr. G. T. Kohman of the Bell Telephone Laboratories has prepar: 

list of such substances that can be mixed, any high one with any low o1 


A 


wom © 





Fig. 2.39-—The refractos« ope 


obtain a fluid satisfactory in all these respects. The following mixtures 
substances are taken from his data. 


| 
Refrac- | Mixture | Flash 
tive Parts | Density P 


Supplier 
Index | by wt 


Point } 
Dimethy! phthalate , ' ‘ta 193 . | Monsanto Chem. Co 
« monochlor naphthalene 3 a 194 oe Bakelite Co. 


Dimethy] phthalate 73 | 1.193 »sop | Monsanto Chem. Co 
Dichlor naphthalene (solid 3+ oo iL Hooker Chem. Co. 
at toom temp.) 


Decalin | 3 | | seno_, | Dupont Co. 
Dowtherm : Dow Chem. Co. 


Kerosene 
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\n immersion fluid for Rochelle salt can be made by mixing decalin with 
inv of the other substances. For example 6 mixture of 34 parts of dimethy1 


phthalate and 66 parts of decalin should give the necessary index 1.495. 


210 APPENDIX 
THE RING ECCENTRICITY CORRECTION 


1 angle a, from the optic axis to- 


Referring to Fig. 2.41, we see that, at a 


vards the plate thickness direction the phase relation is, by Equations (2 
ind (2.9): 


00917 ¢’ sin” ay 





1 = 3 
A cos(6 — Q,) 
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Fig. 2.40—Refraction correction nomograph 
while at an angle a2 away, it is: 
, 00917 t' sin’ ae 
No = - 
A cos(6 + ae) 
Whence, if these are equal, we have: 
sin ay sin’ ae , 
= (2.11) 


cos(é — a) cos(6 + a) 
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These points, in the conoscope field, being of equal phase are parts 
same ring, and if matched to a pair of reticule lines, the optic axis is off + heir 
center line by an angle e, where 


If the separation of the reticule lines corresponds to an angle 2d, we se 


that 
_ at a 











RETICULE 











ill 


Fig. 2.41—Off center cross hairs 


so that 
a =dt+e and a=d-—e 
The reading of the dial will be, at this match 
R= 36 +a+8- a) =58-7 >" 


Hence we have 


sin(d +e) _ 4/28 m @ 


sin(d — e) cos(R + d) 
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ve Ot a For a given value of d, we can plot the values of e as a function of 6 as given 

Off their by (2.12). This plot is a chart of corrections to be added to the readings R 
to find the true angle 6. 

Examination of Equation (2.12) shows us that the correction e is indepen- 

dent of the thickness ¢’ and even of the birefringence; hence, the chart 

Pia could serve for all uniaxial crystals. Equation (2.12) can be given an ap- 


proximate solution: 


1820 tan? d tan R minutes 


II 


é 


\| 


= 30.3 tan? d tan R degrees (2.13) 











Fig. 2.42—Refraction at a surface 


For R not more than 60 degrees and d not more than 15 degrees the error 
in é is not more than 5 minutes. Figure 2.33 is a chart of these corrections 
computed from the more exact equation (2.12). 


ANNULLING THE RING ECCENTRICITY CORRECTION BY MEANS OF THE 
REFRACTION CORRECTION 


lhe difference in quartz index and liquid index is 


and by the law of refraction: 


Nzesin R = Ng sin (R+e) 








262 BELL SYSTEM TECHNICAL JOURNAL 


whence 


sin € 
tan R 


—-H=WN, — (1 — cose) 
and if € is small 
N, sine 


Y= 
tan R 


From Eq. 2.13, e = 0.530 tan? d tan R radians, and putting this in the 
tion for #7 we find that the correction for ring eccentric ity approxin 


annuls the correction for refraction if 


HT = — .530 tan? d 


For d = 10° this gives // - 0255, that is, a fluid index of 1.5207. 














A Note on the Transmission Line Equation 
in Terms of Impedance 


By J. R. PIERCE 


NCREASED familiarity derived in handling Maxwell’s equations, 
] especially in connection with problems arising at very high frequencies, 
has resulted in a variety of forms for expressing certain laws and behavior. 
Especially, work by Schelkunoff in extending the impedance concept! shows 
that impedance can be quite as general and exact a means for expressing 
electromagnetic relations as are current, voltage, electric and magneti 
fields, and vector and scalar potentials. 

In reformulating certain problems in terms of impedance the content and 
iltimate solution must of course be equivalent. There may, however, be 
a considerable change of procedure and sometimes a simplification. For 
nstance, in many cases a single impedance condition can replace the usual 
two boundary conditions for voltage and current. 

One very simple case in which it is perhaps easiest to deal directly with 
impedance is in the derivation of the transmission line equation on a dis- 
tributed constant basis. In the usual derivation, two linear second order 
differential equations are obtained, one for voltage and one for current. 
The impedance, in terms of which the engineer expresses many of his results, 
is obtained as a ratio from solutions for voltage and current. In treating 
the transmission line from the impedance point of view, without dealing 
with currents and voltages, a first order non-linear differential equation in 
terms of impedance and distance is obtained. This impedance equation 
is a Ricatti equation and could be obtained from the usual line equations. 
It is simpler, however, to derive it directly. 

As the principal interest of such a treatment lies in the method and in 
the fact that the line may be tapered, rather than in losses, the derivations 
will be carried out for lossless lines. Losses can be taken into account by 
allowing the inductance per unit length, 1, and the capacitance per unit 
length, C, to become complex quantities. 

Consider the section of line dy long, shown in the figure, having an 
inductance L dx and a capacitance C dx. We can write immediately 


bx + dZ = Zs+d2 
1 
jwC dx + 1/Z, 
= Z, + jo[L — CZ%| dx. 1) 


= jwL dx + 


'“The Impedance Concept and Its Application to Problems of Reflection, Refraction, 
Shielding, and Power Absorption,” B.S.7.J. Vol. 17, pp. 17-48, January, 1938 
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Dropping the subscript «, the differential equation for the line in te 
the impedance Z may be written* 


R 


R = (L/C)'" 
= (icy 


R is the nominal characteristic impedance, and v is the nominal phas 
velocity, which is constant for many tapered lines with the same dic|ectri 
material separating the conductors throughout their length. In such lines, 
if the dielectric is air or vacuum, 2 is c, the velocity of light. 

It should not be surprising that (2) is of the first order. Although there 
are two boundary conditions, the impedances terminating the right and 
left ends of the line, there are two impedances, that looking toward the right 
and that looking toward the left. The impedance looking toward the right 


Lax 





Fig. 1 


is unaffected by the left end termination, and that looking toward the left 
is unaffected by the right end termination. 

As R is real, it may be seen from (2) that the only case in which the im 
pedance can equal the nominal characteristic impedance R at all points 
is for R constant. This tells us that the characteristic impedance of an) 
lossless tapered line is complex. For very gradually tapering lines th 
characteristic impedance differs from the nominal characteristic impedanc 
principally by a small imaginary component. 

The simplest solution of (2) is of course that for a uniform line, with R 
a constant which will be called R,. In this case (2) can be integrated 
directly, giving the familiar result 


a tanh (jwx/v + K) 
Ro 


2 It is interesting to note that the equation for admittance Y can be obtained by r 
placing Z by Y and R by (1/R) = G in (2). 
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t . / : 
Dr. L. A. MacColl has pointed out to the writer that (2) is the same as 
the electrostatic electron optical equation for paraxial rays. To reduce (2 
to the standard form: 
) 
jw dx d 
= = dz ¢ 
3 Rv (6) 
(4 —R’ = H(z) (7) 
dZ 2 
| phas = H(z) +2 (8) 
icles r dz 
h lines [he electron optical equation for paraxial rays is 
) rr 2 
h there dl = 3 Ka + ‘ig (9) 
At and dz 16 L V(z) 
he right V'(z 
» right r=C- ( (10) 
1€ righ 4V(z) 


Here z is distance along the axis, V(s) is potential on the axis, and C is 
convergence, or the inverse of the focal distance. 

It would seem, then, that from each solution of an electron optical 
problem, a solution of a tapered line problem could be found, and vice versa. 

While it cannot be claimed that anything new has entered the transmis- 
sion line equation in expressing it in terms of impedance, it does seem that 
the approach may be stimulating in uncovering hitherto neglected material 
and analogies. 
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Abstracts of Technical Articles by Bell System Authors 


Electronics in Telephony.) Frank A. Cowan. The historical dev: 
ment of the use of electronic devices by the telephone system is review 
showing how long distance telephony has grown with the increased us 
and improvements in, electronics. The number of telephone repeaters 
grown from 16 in 1908 to 123,000 in 1942 and carrier circuit mileage 
grown from 2,000 in 1920 to 2,300,000 in 1942, while copper usage per ci! 
mile has decreased from 400 pounds in 1910 to less than 70 in 1942. 

A transcontinental telephone connection has grown from an open 
circuit with a total loss, less repeaters, of less than 75 db (1915) to a pres 
day cable circuit operating at carrier frequencies, which may have a t 
loss of over 10,000 db. The problem of matching enormous amplificat 
to compensate for huge losses with a precision of one or two db wa 
difficult one, which was solved by electronic techniques. The ampliii 
tion necessary to compensate for the high losses on the cable layout 
entail some 200 repeaters utilizing a total of more than 600 vacuum tul: 
in tandem. The automatic regulation and control of the amplificatio: 
accomplished by electronic devices, giving to the present day circuits 
stability unequalled in the days before the vacuum tube. 

There is available, except for the War, radiotelephone service to » 
foreign countries and overseas areas, and ocean liners at sea, and to boat 
in coastal and inland waters. 

Such widespread dependence on vacuum tubes has stimulated resear 
and design to achieve long life and a high degree of uniformity, stabi 
and reliability. Among interesting future possibilities are transocea! 
cables, the use of higher frequencies providing broader bands and large: 
numbers of circuits over a given path, and further application of remote ai 
unattended stations. 


Deionization Considerations in a Harmonic Generator Employing a Ga 
Tube Switch? WittiAM G. SHEPHERD. A description is given ol 
experimental investigation of the properties of a thyratron operating as 
high-frequency switch in a circuit which permitted the generation o 
wide band of harmonics. ‘The experiments indicate that there is an operat 
ing frequency below which no difficulties in deionization occur and ab: 

' Electronics, March 1943. 

> Proc. 1.R.E., February 1943. 
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ABSTRACTS OF TECHNICAL ARTICLES 


h stable operation requires that the grid potential fulfill certain cot 
‘ions dependent upon the frequency, wave form of r 
ut constants. It has been found possible to operate certain 
vratrons at switching frequencies as high as several h 
el second. For these higher frequent ies the deionization ol the bes 1s 
devi complete but normal switching behavior is obtained 
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Insas For the purposes of record and assistance to librarians, 
hone and for the information of subscribers, it is to be noted 
IS Ol that there was no April 1943 issue of the Bell System 
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